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Abstract 

A general geometrical framework of nonequilibrium thermodynamics is devel- 
oped. The notion of macroscopically definable ensembles is developed. The thesis 
about macroscopically definable ensembles is suggested. This thesis should play the 
same role in the nonequilibrium thermodynamics, as the Church- Turing thesis in 
the theory of computability. The primitive macroscopically definable ensembles are 
described. These are ensembles with macroscopically prepared initial states. The 
method for computing trajectories of primitive macroscopically definable nonequi- 
librium ensembles is elaborated. These trajectories are represented as sequences of 
deformed equilibrium ensembles and simple quadratic models between them. The 
primitive macroscopically definable ensembles form the manifold in the space of en- 
sembles. We call this manifold the film of nonequilibrium states. The equation for 
the film and the equation for the ensemble motion on the film are written down. 
The notion of the invariant film of non-equilibrium states, and the method of its 
approximate construction transform the the problem of nonequilibrium kinetics into 
a series of problems of equilibrium statistical physics. The developed methods al- 
low us to solve the problem of macro-kinetics even when there are no autonomous 
equations of macro-kinetics. 
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Introduction 



The goal of this paper is to discuss the nonhnear problems of irreversibility, and to revise 
the previous attempts to solve them. The interest to the problem of irreversibility persists 
during decades. It has been intensively discussed in the past, and nice accounts of these 
discussions can be found in the literature (see, for example, PQ 121111 El)- Here, we intend 
to develop a more geometrical viewpoint on the subject. The paper consists of two 
parts. First, in section^ we discuss in an informal way the origin of the problem, and 
demonstrate how the basic constructions arise. Second, in section |21 we give a consistent 
geometric formalization of these constructions. Our presentation is based on the notion of 
the natural projection introduced in sectional We discuss in detail the method of natural 
projector as the consistent formalization of Ehrenfest's ideas of coarse-graining. 

In section 13.21 we introduce a one-dimensional model of nonequilibrium states. In the 
background of many derivations of nonequilibrium kinetic equations there is present the 
following picture: Above each point of the quasiequilibrium manifold there is located a 
huge subspace of nonequilibrium distributions with the same values of the macroscopic 
variables, as in the quasiequilibrium point. It is as if the motion decomposes into two 
projections, above the point on the quasiequilibrium manifold, and in the projection on 
this manifold. The motion in each layer above the points is highly complicated, but fast, 
and everything quickly settles in this fast motion. 

However, upon a more careful looking into the motions which start from the quasiequi- 
librium points, we will observe that, above each point of the quasiequilibrium manifold 
it is located just a single curve, and all the nonequilibrium (not-quasiequilibrium) states 
which come into the game form just a one-dimensional manifold. 

The novel approach developed in section |21 allows to go beyond limitations of the short 
memory approximations through a study of stability of the quasiequilibrium manifold. 

The one-dimensional models of nonequilibrium states form a film of nonequilibrium 
states over the quasiequilibrium manifold. In section H) we present a collection of methods 
for the film construction. One of the benefits from this new technic is the possibility to 
solve the problem of macro-kinetic in cases when there are no autonomous equations of 
macro-kinetic for moment variables. The notion of the invariant film of non-equilibrium 
states, and the method of its approximate construction transform the the problem of 
nonequilibrium kinetics into a series of problems of equilibrium statistical physics. To 
describe a dynamics of nonequilibrium ensemble one should find series of deformed equi- 
librium ensembles. 

In Appendix a short presentation of the method of invariant manifold for kinetic 
problems is given. 

The most important results of the paper are: 

1. The notion of macro scopically definable ensembles is developed. 

2. The primitive macro scopically definable ensembles are described. 
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3. The method for computing trajectories of primitive macroscopically definable nonequi- 
hbrium ensembles is elaborated. These trajectories are represented as series of de- 
formed equilibrium ensembles and simple quadratic models between them. 

Let us give here an introductory description of these results. 

The notion of macroscopically definable ensembles consists of three ingredients: 

1. Macroscopic variables, the variables which values can be controlled by us; 

2. Quasiequilibrium state, the conditional equilibrium state for fixed values of the 
macroscopic variables; 

3. Natural dynamics of the system. 

We use the simplest representation of the control: For some moments of time we 
fix some values of the macroscopic variables (of all these variables, or of part of them; 
for the whole system, or for macroscopically defined part of it; the current "natural", or 
some arbitrary value of these variables), and the system obtains corresponding conditional 
equilibrium state. We can also keep fixed values of some macroscopic variables during a 
time interval. 

These control operations are discrete in time. The continuous control can be cre- 
ated by the closure: the limit of a sequence of macroscopically definable ensembles is 
macroscopically definable too. 

The role of the macroscopic variables for the irreversibility problem became clear to M. 
Leontovich and J. Lebowitz several decades ago jSl El 13 E] • But this was the first step. 
Now we do need the elaborate notion of ensembles which can be obtained by macroscopic 
tools. The Maxwell Demon gives the first clear picture of a difference between macroscopic 
and microscopic tools for ensembles control (there are books devoted to analysis of this 
Demon piU ITT]). Nevertheless, the further step to the analysis of the notion of macroscopic 
definability in context of constructive transition from microdynamics to macrokinetics 
equations had not done before the paper [12j. Our analysis pretends to be an analogue 
of the Church- Turing thesis ^3 E| • This thesis concerns the notion of an effective (or 
mechanical) method in mathematics. As a "working hypothesis". Church proposed: A 
function of positive integers is effectively calculable only if recursive. 

We introduce a class of "macroscopically definable ensembles" and formulate the the- 
sis: An ensemble can be macroscopically obtained only if macroscopically definable in 
according to the introduced notion. This is a thesis about success of the formalization, as 
the Church- Turing thesis, and nobody can prove or disprove it in rigorous sense, as well, 
as this famous thesis. 

Another important new notion is the "macroscopically definable transformation" of 
the ensemble: If one get an ensemble, how can he transform it? First, it is possible just 
let them evolve, second, it can be controlled by the macroscopic tools on the defined way 
(it is necessary just to keep values of some macroscopic variables during some time). 
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The primitive macro scopically definable ensembles aiQ ensembles with quasiequihbrium 
initial states and without further macroscopic control. These ensembles are prepared 
macroscopically, and evolve due to natural dynamics. The significance of this class of 
ensembles is determined by the hypothesis about the primitive macroscopically definable 
ensembles: Any macroscopically definable ensemble can be approximated by primitive 
macroscopically definable ensembles with appropriate accuracy. Now we have no other 
effective way to decribe the nonequilibrium state. 

The primitive macroscopically definable ensembles form the manifold in the space of 
ensembles. We call this manifold the "film of nonequilibrium states". The equation for 
the film and the equation for the ensemble motion on the film are written down. 

The film of nonequilibrium states is the trajectory of the manifold of initial quasiequi- 
hbrium states due to the natural (microscopic) dynamics. For every value of macroscopic 
variables this film gives us a curve. The curvature of this curve defines kinetic coefficients 
and entropy production. 

The main technical problem is the computation of this curve for arbitrary value of the 
macroscopic variables. We represent it as a sequence of special points and second-order 
polynomial (Kepler) models for trajectory between these points. The method elaborated 
for the computation is the further development of the method for initial layer problem in 
the Boltzmann kinetics P^ ll6j. For dissipative Boltzmann microkinetics it was sufficient 
to use the first-order models (with or without smoothing). For conservative microkinetics 
it is necessary to use the highest-order models. Application of this method to the lattice 
kinetic equations gave the following possibilities: 

• To create the Lattice-Boltzmann kinetics with if-theorem [T7] : 

• To transform the Lattice-Boltzmann method into the numerically stable computa- 
tional tool for fiuid fiows and other dissipative systems out of equilibrium jl8j : 

• To develop the Entropic Lattice Boltzmann method as a starting basis for the for- 
mulation of a new class of turbulence models based on genuinely kinetic principles 

m- 

Now we extend the method elaborated for dissipative systems HE] to the higher- 
order models for conservative systems. The iteration method for improvement of obtained 
approximations is proposed too. It is a version of the Method of invariant manifold for 
kinetic problems, developed in the series of papers [201 1^ (the almost exhaustive review 
of these works can be find in the paper (221) • summary of this method is given in 
Appendix. 

The constructing of the method of physically consistent computation is the central 
part of our paper. It is neither a philosophical opus, nor only discussion of foundations 
of science. 

The main results of this paper were presented in the talk given on the First Mexican 
Meeting on Mathematical and Experimental Physics, Mexico City, September 10-14, 2001, 
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and in the lectures given on the V Russian National Seminar "Modeling of Nonequilibrium 
systems", Krasnoyarsk, October 18-20, 2002 gHl. 

1 The problem of irreversibility 

1.1 The phenomenon of the macroscopic irreversibihty 

The "stairs of reduction" (Fig. lead from the reversible microdynamics to irreversible 
macrokinetics. The most mysterious is the first step: the appearance of irreversibility. 

The best way to demonstrate the problem of irreversibility is the following Gedanken- 
experiment Let us watch the movie: It's raining, people are running, cars rolling. Let 
us now wind this movie in the opposite direction, and we will see a strange and funny 
picture: Drops of the rain are raising up to the clouds, which next condensate into the 
vapor on the pools, on the surfaces of rivers, people run with their backs forward, cars 
behave also quite strange, and so forth. This cannot be, and we "know" this for sure, we 
have never seen anything like this in our life. Let us now imagine that we watch the same 
movie with a mag nitude of 10*^ - 10^ so that we can resolve individual particles. And all 
of the sudden we discover that we cannot see any substantial difference between the direct 
and the reverse demonstration: Everywhere the particles are moving, colliding, reacting 
according to the laws of physics, and nowhere there is a violation of anything. We cannot 
tell the direct progressing of the time from the reversed. So, we have the irreversibility of 
the macroscopic picture under the reversibility of the microscopic one. 

Rain, people, cars - this all is too complicated. One of the most simple examples of 
the irreversible macroscopic picture under the apparent reversibility of the microscopic 
picture (the "thermal ratchet") is given by R. Feynman in his lectures on the character 
of physical law f23j. We easily label it as self-evident the fact that particles of different 
colors mix together, and we would see it as a wonder the reverse picture of a spontaneous 
decomposition of their mixture. However, itself an appreciation of one picture as usual, 
and of the other as unusual and wonderful - this is not yet the physics. It is desirable to 
measure somehow this transition from order to disorder. 

1.2 Phase volume and dynamics of ensembles 

Let there be n blue and n white particles in a box, and let the box is separated in two 
halves, the left and the right. Location of all the particles in the box is described by the 
assembly of 2n vectors of locations of individual particles. The set of all the assemblies 
is a "box" in the 6n-dimensional space. A point in this 6n- dimensional box describes a 
configuration. The motion of this point is defined by equations of mechanics. 

"Order" is the configuration in which the blue particles are all in the right half, and all 
the white particles are in the left half. The set of all such configurations has a rather small 
volume. It makes only (1/2)^" of the total volume of the Gn-dimensional box. If n = 10, 
this is of the order of one per million of the total volume. It is practically unthinkable 
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to get into such a configuration by chance. It is also highly improbable that, by forming 
more or less voluntary the initial conditions, we can observe that the system becomes 
ordered by itself. Prom this standpoint, the motion goes from the states of "order" to the 
state of "disorder" , just because there are many more states of "disorder" . 

However, we have defined it in this way. The well known question of where is more 
order, in a fine castle or in a pile of stones, has a profound answer: It depends on which 
pile you mean. If "piles" are thought as all configurations of stones which are not castles, 
then there are many more such piles, and so there is less order in such a pile. However, 
if these are specially and uniquely placed stones (for example, a garden of stones), then 
there is the same amount of order in such a pile as in the fine castle. Not a specific 
configuration is important hut an assembly of configurations embraced by one notion. 

This transition from single configurations to their assemblies (ensembles) play the 
pivotal role in the understanding of irreversibility: The irreversible transition from the 
ordered configuration (blue particles are on the right, white particles are on the left) 
to the disordered one occurs simply because there are many more of the disordered (in 
the sense of the volume). Here, strictly speaking, we have to add also a reference to the 
Liouville theorem: The volume in the phase space which is occupied by the ensemble does 
not change in time as the mechanical system evolves. Because of this fact, the volume V 
is a good measure to compare the assemblies of configurations. However, more often the 
quantity In V is used, this is called the entropy. 

The point which represents the configuration, very rapidly leaves a small neighborhood 
and for a long time (in practice, never) does not come back into it. In this, seemingly 
idyllic picture, there are still two rather dark clouds left. First, the arrow of time has not 
appeared. If we move from the ordered initial state (separated particles) backwards in 
time, then everything will stay the same as when we move forward in time, that is, the 
order will be changing into the disorder. Second, let us wind the film backwards, let us 
shoot the movie about mixing of colored particles, and then let us watch in the reverse 
order their demixing. Then the initial configurations for the reverse motion will only seem 
to be disordered. Their "order" is in the fact that they were obtained from the separated 
mixture by letting the system to evolve for the time t. There are also very few such 
configurations, just the same number as of the ordered (separated particles) states. If we 
start with these configurations, then we obtain the ordered system after the time t. Then 
why this most obvious consequence of the laws of mechanics looks so improbable on the 
screen? Perhaps, it should be accepted that states which are obtained from the ordered 
state by a time shift, and by inversion of particle's velocities (in order to initialize the 
reverse motion in time), cannot be prepared by using macroscopic means of preparation. 
In order to prepare such states, one would have to employ an army of Maxwell's daemons 
which would invert individual velocities with sufficient accuracy (here, it is much more 
into the phrase "sufficient accuracy" but this has to be discussed separately and next 
time) . 

Por this reason, we lump the distinguished initial conditions, for which the mixture 
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decomposes spontaneously ("piles" of special form, or "gardens of stones") together with 
other configurations into macro scopically definable ensembles. And already for those en- 
sembles the spontaneous demixing becomes improbable. This way we come to a new 
viewpoint: (i). We cannot prepare individual systems but only representatives of ensem- 
bles, (ii) We cannot prepare ensembles at our will but only "macroscopically definable 
ensembles" . What are these macroscopically definable ensembles? It seems that one has 
to give some constructions, the universality of which can only be proven by the time and 
experience. 

1.3 Macroscopically definable ensembles and quasiequilibria 

The main tool in the study of macroscopically definable ensembles is the notion of the 
macroscopic variables, and of the quasiequilibria. In the dynamics of the ensembles, the 
macroscopic variables are defined as linear functionals (moments) of the density distri- 
bution of the ensemble. Macroscopic variables M usually include hydrodynamic fields, 
density of particles, densities of momentum, and density of the energy, also the list may 
include stress tensor, reaction rates and other quantities. In the present context, it is 
solely important that the list the macroscopic variables is identified for the system under 
consideration. 

A single system is characterized by a single point x in the phase space. The ensemble 
of the systems is defined by the probability density F on the phase space. Density F must 
satisfy a set of restrictions, the most important of which are: Nonnegativity, F[x) > 0, 
normalization, 

[ F{x)dV{x) = 1, (1) 
Jx 

and that the entropy is defined, that is, there exists the integral, 

S{F) = - I F{x)hiF{x)dV{x). (2) 
Jx 

(Function FlnF is continuously extended to zero values of F: OlnO = 0). Here dV{x) is 
the invariant measure (phase volume). 

The quasiequilibrium ensemble describes the "equilibrium under restrictions". It is 
assumed that some external forcing keeps the given values of the macroscopic variables 
M, with this, "all the rest" comes the corresponding (generalized) canonic ensemble F 
which is the solution to the problem: 

S{F) max, M(F) = M. (3) 

where S{F) is the entropy, M{F) is the set of macroscopic variables. 

The thesis about the macroscopically definable ensembles. Macroscopically 
definable ensembles are obtained as the result of two operations: 

(i). Bringing the system into the quasiequilibrium state corresponding to either the 
whole set of the macroscopic variables M, or to its subset. 
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(ii). Changing the ensemble according to the microscopic dynamics (due to the Liou- 
ville equation) during some time t. 

These operations can be apphed in the interchanging order any number of times, and 
for arbitrary time segments t. The hmit of macroscopically definable ensembles will also 
be termed the macroscopically definable. One always starts with the operation (i). 

In order to work out the notion of macroscopic definability, one has to pay more 
attention to partitioning the system into subsystems. This involves a partition of the phase 
space X with the measure dV on it into a direct product of spaces, X = Xi x X2 with 
the measure dVidV2. To each admissible ("macroscopic") partition into sub-systems, it 
corresponds the operation of taking a "partial quasiequilibrium" , applied to some density 
Fo(xi,X2): 

S{F) max, (4) 
M{F) = M, I F{xi,X2)dV2{x2)= f Fo{xi,X2)dV2{x2). 

where M is some subset of macroscopic variables (not necessarily the whole list of the 
macroscopic variables). In Eq. (jD), the state of the first subsystem is not changing, 
whereas the second subsystem is brought into the quasiequilibrium. In fact, the problem 
(0)) is a version of the problem with additional "macroscopic variables" , 

/ F{xuX2)dV2{x2). (5) 

JX2 

The extended thesis about macroscopically definable ensembles allows to use also 
operations with only one restriction: The initial state should be the "true quasiequi- 
librium" that is, macroscopic variables related to all possible partitions into subsystems 
should appear only after the sequence of operations has started with the solution to the 
problem Q for some initial M. This does not exclude a possibility of including oper- 
ators (0) into the list of the basic macroscopic variables M. The standard example of 
such an inclusion are few-body distribution functions treated as macroscopic variables in 
derivations of kinetic equations from the Liouville equation. 

Irreversibility is related to the choice of the initial conditions. The extended set of 
macroscopically definable ensembles is thus given by three objects: 

(i) . The set of macroscopic variables M which are linear (and, in an appropriate 
topology, continuous) mappings of the space of distributions onto the space of values of 
the macroscopic variables; 

(ii) . Macroscopically admissible partitions of the system into sub-systems; 

(iii) . Equations of microscopic dynamics (the Liouville equation, for example). 

The choice of the macroscopic variables and of the macroscopically admissible par- 
titions is a distinguished topic. The main question is: what variables are under the 
macroscopic control? Here the macroscopic variables are represented as formal elements 
of the construction, and the arbitrariness is removed only at solving specific problems. 
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Usual we can postulate some properties of macroscopic variables, for example, symmetry 
with respect to any permutation of equal micro-particles. 

We have discussed the prepared ensembles. But there is another statement of problem 
too: Let us get an ensemble. The way how we get it may be different and unknown, for 
example, some demon or oracle^ can give it to us. How can we transform this ensemble by 
macroscopic tools? First, it is possible just let them evolve, second, it can be controlled 
by the macroscopic tools on the defined way (it is necessary just to keep values of some 
macroscopic variables during some time). 

The thesis about the macroscopically definable transformation of ensem- 
bles. Macroscopically definable transformation of ensembles are obtained as the result of 
two operations: 

(i) . Bringing the system into the quasiequilibrium state corresponding to either the 
whole set of the macroscopic variables M, or to its subset. 

(ii) . Changing the ensemble according to the microscopic dynamics (due to the Liou- 
ville equation, for example) during some time t. 

These operations can be applied in the interchanging order any number of times, and 
for arbitrary time segments t. The limit of macroscopically definable transformations 
will also be termed the macroscopically definable. The main difference of this definition 
(macroscopically definable transformation) from the definition of the macroscopically de- 
finable ensembles is the absence of restriction on the initial state, one can start from 
arbitrary ensemble. 

The class of macroscopically definable ensembles includes the more simple, but im- 
portant class. Let us reduce the macroscopic control to preparation of quasiequilibrium 
ensemble: we just prepare the ensemble by macroscopic tools and then let them evolve 
due to natural dynamics (Liouville equation, for example). Let us call this class the prim- 
itive macroscopically definable ensembles. These ensembles appear as results (for t > 0) 
of motions which start from the quasiequilibrium state (at t = 0). The main technical 
results of our work concern the computation of the manifold of primitive macroscopically 
definable ensembles for a given system. 

The importance of this class of ensembles is determined by the following hypothesis: 
The hypothesis about the primitive macroscopically definable ensembles. Any 
macroscopically definable ensemble can be approximated by primitive macroscopically 
definable ensembles with appropriate accuracy. In some limits we can attempt to say: 
"with any accuracy" . Moreover, this hypothesis with "arbitrary small accuracy" can be 
found as the basic but implicit foundation of all nonequilibrium kinetics theories which 
pretend to demonstrate a way from microdymamics to macrokinetics, for example in 
Zubarev nonequilibrium statistical operator theory 3j, etc. This hypothesis allow to 
describe nonequilibrium state as a result of evolution of quasiequilibrium state in time. 
Now we have no other way to decribe the nonequilibrium state^ 

^In the theory of computation, if there is a device which could answer questions beyond those that a 
Turing machine can answer, then it would be called oracle. 

^ There exists a series of papers with discussion of Hamiltonian systems in so-called force thermostat, 
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The hypothesis about the primitive macroscopically definable ensembles is real hy- 
pothesis, it can hold for different systems with different accuracy, it can be true or false. 
In some limits the set of primitive macroscopically definable ensembles can be dense in 
the set of all macroscopically definable ensembles, or can be not dense, etc. There is 
the significant difference between this hypothesis and the thesis about macroscopically 
definable ensembles. The thesis can be accepted, or not, the reasons for acceptance can 
be discussed, but nobody can prove or disprove the definition, even the definition of the 
macroscopically definable ensembles. 

1.4 Irreversibility and initial conditions 

The choice of the initial state of the ensemble plays the crucial role in the thesis about the 
macroscopically definable ensembles. The initial state is always taken as the quasiequi- 
librium distribution which realizes the maximum of the entropy for given values of the 
macroscopic variables. The choice of the initial state splits the time axis into two semi- 
axes: moving forward in time, and moving backward in time, the observed non-order 
increases (the simplest example is the mixing of the particles of different colors). 

In some works, in order to achieve the "true nonequilibrium" , that is, the irreversible 
motion along the whole time axis, the quasiequilibrium initial condition is shifted into 
— oo in time. This trick, however, casts some doubts, the major being this: Most of 
the known equations of macroscopic dynamics which describe irreversible processes have 
solutions which can be extended backwards in time only for finite times (or cannot be 
extended at all). Such equations as the Boltzmann kinetic equation, diffusion equation, 
equations of chemical kinetics and like do not allow for almost all their solutions to be 
extended backward in time for indefinitely long. All motions have a "beginning" beyond 
which some physical properties of a solution will be lost (often, positivity of distributions), 
although formally solutions may even exist, as in the case of chemical kinetics. 

1.5 Weak and strong tendency to equilibrium, shaking and short 
memory 

One aspect of irreversibility is the special choice of initial conditions. Roughly speaking, 
the arrow of time is defined by the fact that the quasiequilibrium initial condition was in 
the past. 

This remarkably simple observation does not, however, exhaust the problem of transi- 
tion from the reversible equations to irreversible macroscopic equations. One more aspect 

or, in particular, isokinetic thermostat (see, for example, the review of D. Ruelle |21j). These thermostats 
were invented in computational molecular dynamics for acceleration of computations, as a technical trick. 
From physical point of view this is a theory about a friction of particles on the space, the "ether friction" 
(the "ether" is a theoretical substrate in the ancient physics). Of course, this theory is mathematically 
consistent and perhaps it may be useful as the theory of special computations methods, but a bridge 
between this theory and physics is desirable. 
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deserves a serious consideration. Indeed, distribution functions tend to the equilibrium 
state according to macroscopic equations in a strong sense: deviations from the equihb- 
rium tends to zero in the sense of most relevant norms (in the sense, for example, 
or even uniformly). On the contrast, for the Liouville equation, tendency to equilibrium 
ocures (if at all) only in the weak sense: mean values of sufficiently "regular" functions 
on the phase space do tend to their equilibrium values but the distribution function itself 
does not tend to the equilibrium with respect to any norm, not even point-wise. This 
is especially easy to appreciate if the initial state has been taken as the equipartition 
over some small bounded subset of the phase space (the "phase drop" with small, but 
non-zero volume). This phase drop can mix over the phase space, but for all the times 
it will remain "the drop of oil in the water", the density will always be taking only two 
values, and p > 0, and the volume of the set where the density is larger than zero 
will not be changing in time, of course. So, how to arrive from the weak convergence 
(in the sense of the convergence of the mean values), to the strong convergence (to the 

or to the uniform convergence, for example)? In order to do this, there are two basic 
constructions: The coarse-graining (shaking) in the sense of Ehrenfests', and the short 
memory approximation. 

The idea of coarse-graining dates back to P. and T. Ehrenfests, and it has been most 
clearly expressed in their famous paper of 1911 Ehrenfests considered a partition of 
the phase space into small cells, and they have suggested to alter the motions of the phase 
space ensemble due to the Liouville equation with "shaking" - averaging of the density 
of the ensemble over the phase cells. In the result of this process, the convergence to the 
equilibrium becomes strong out of the weak. It is not difficult to recognize that ensembles 
with constant densities over the phase cells are quasiequilibria; corresponding macroscopic 
variables are integrals of the density over the phase cells ("occupation numbers" of the 
cells). This generalizes to the following: alternations of the motion of the phase ensemble 
due to microscopic equations with returns to the quasiequilibrium manifold, preserving 
the values of the macroscopic variables. It is precisely this construction which serves for 
the point of departure for many of the constructions below. 

Another construction is the short memory approximation. The essence of it is the 
following: If one excludes microscopic variables and assumes quasiequilibrium initial con- 
ditions, then it is possible to derive integro-differential equations with retardation for 
the macroscopic variables (the way to do this is not unique). The form of the resulting 
equations is approximately this: 



where K{t,t') is an operator (generally speaking, nonlinear) acting on M{t'). Once this 
equation is obtained, one assumes that the kernels of these integro-differential equations 
decay at a sufficiently high rate into the past (for example, exponentially, as ||-ft'(t, t') [M{t')] \\ < 
exp{ — (t — t')/r}||M(t')||). This can be interpreted in the spirit of Ehrenfests': Every mo- 
tion which has begun sufficiently recently (the "memory time" r before now) can be 
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regarded as being started from the quasiequilibrium. Thus, after each time r has elapsed, 
the system can be shaken in the sense of Ehrenfests - the result should not differ much. 

1.6 The essence of irreversibility in two words 

(i) The direction of the arrow of time is defined by the fact that only "macroscopically 
definable ensembles" can be taken as initial conditions, that is, only quasiequilibrium 
ensembles and what can be obtained from them when they are exposed to the true mi- 
croscopic dynamics, or when partial quasiequilibria are taken in positive time. We are 
created in such a way that we prepare and control (in part) the present, and observe what 
happens in the future. (In a sense, this is a definition of the subjective time). 

(ii) Microscopic dynamics can give only the weak convergence to the equilibrium, conver- 
gence of mean values. Macroscopic variables tend to the equilibrium in the strong sense. 
The passage from micro to macro occurs here with the help of Ehrenfests' coarse-graining 
procedure or its analogs. 

One might feel uneasy about the second of these points because the procedure of 
coarse-graining is not the result of the equations of motion, and therefore it is somehow 
voluntary. The only hope to lift this arbitrariness is that it may well happen that, in the 
limit of a very large number of particles, the perturbation caused by the coarse-graining 
can be made arbitrary small, for example, by increasing the time interval between coarse- 
graining. 

1.7 Equivalence between trajectories and ensembles in the ther- 
modynamic limit 

In the preceding sections we were speaking about the dynamics of ensembles. This ap- 
parently contradicts the fact that the dynamics of a classical system goes along a single 
trajectory. Two arguments make it possible to proceed from the trajectories to ensembles: 

(i) High sensitivity of trajectories to external perturbations when the number of par- 
ticles is large. Arbitrary weak noise results in the stochastization of the motion. 

(ii) In the thermodynamic limit, it is possible to partition the system into an arbitrary 
large number of small but still macroscopic sub-systems. Initial conditions in the sub- 
systems are independent from one sub-system to another, and they cannot be assigned 
completely voluntary but are taken from some distribution with a fixed sum of mean values 
(an analog of the macroscopic definability of ensembles). For spatially inhomogeneous 
systems, such small but still macroscopic subsystems are defined in small and "almost 
homogeneous" volumes. 

1.8 Subjective time and irreversibility 

In our discussion, the source of the arrow of time is, after all, the asymmetry of the 
subjective time of the experimentalist. We prepare initial conditions, and after that we 
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watch what will happen in the future but not what happened in the past. Thus, we 
obtain kinetic equations for specifically prepared systems. How is this related to the 
dynamics of the real world? These equations are applicable to real systems to the extent 
that the reality can be modeled with systems with specifically prepared quasiequilibrium 
initial conditions. This is anyway less demanding than the condition of quasi-staticity 
of processes in classical thermodynamics. For this reason, versions of nonequilibrium 
thermodynamics and kinetics based on this understanding of irreversibility allowed to 
include such a variety of situations, and besides that, they include all classical equations 
of nonequilibrium thermodynamics and kinetics. 

2 Geometrization of irreversibility 
2.1 Quasiequilibrium manifold 

Let E he a linear space, and U G E he a convex subset, with a nonempty interior intU. 
Let a twice differentiable concave functional S he defined in intf/, and let S is continuous 
on U. According to the familiar interpretation, S is the entropy, E is an appropriate space 
of distributions, U is the cone of nonnegative distributions from E. Space E is chosen in 
such a way that the entropy is well defined on U. 

Let L be a closed linear subspace of space E, and m : E E/L he the natural 
projection on the factor-space. The factor-space E/L will further play the role of the 
space of macroscopic variables (in examples, the space of moments of the distribution). 

For each M G intU/L we define the quasiequilibrium, flj- G int?7, as the solution to 
the problem, 

S{f) ^ max, m(/) = M. (6) 

We assume that, for each M G intU/L, there exists the (unique) solution to the problem 
(P). This solution, fl^, is called the quasiequilibrium, corresponding to the value M 
of the macroscopic variables. The set of quasiequilibria fl^ forms a manifold in intf/, 
parameterized by the values of the macroscopic variables M G intU/L (Fig. 

Let us specify some notations: E^ is the adjoint to the E space. Adjoint spaces 
and operators will be indicated by , whereas notation * is earmarked for equilibria and 
quasiequilibria. 

Furthermore, [l,x] is the result of application of the functional / G E'^ to the vector 
X G -E. We recall that, for an operator A : Ei ^ E2, the adjoint operator, : Ef E2 
is defined by the following relation: For any / G E2 and x G Ei, 

[I, Ax] = [A^l,x]. 

Next, DfS{f) G E'^ is the differential of the functional S{f), D^S{f) is the second 
differential of the functional S{f). Corresponding quadratic functional D'^S{f){x,x) on 
E is defined by the Taylor formula. 
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QE Manifold 




Figure 2: Relations between a microscopic state /, the corresponding macroscopic state 
M = m{f), and quasiequilibria fl,j. 
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S{f + x) = S{f) + [DfSif),x] + ^DjS{f){x, x) + o(||xf ). (7) 

We keep the same notation for the corresponding symmetric bihnear form, Dj:S{f){x,y), 
and also for the hnear operator, D'jSlf) : E E^, defined by the formula, 

[DjS{f)x,y] = DjS{f){x,y). 

Here, on the left hand side there is the operator, on the right hand side there is the 
bilinear form. Operator D'jS{f) is symmetric on E, D^S{f)'^ = D'jS{f). 

Concavity of S means that for any x E E the inequality holds, D'jS{f){x,x) < 0; 
in the restriction onto the affine subspace parallel to L we assume the strict concavity, 
DjS{f){x,x) < if a; e L, and X 7^ 0. 

A comment on the degree of rigor is in order: the statements which will be made below 
become theorems or plausible hypotheses in specific situations. Moreover, specialization 
is always done with an account for these statements in such a way as to simplify the 
proofs. 

Let us compute the derivative Dm/m- For this purpose, let us apply the method of 
Lagrange multipliers: There exists such a linear functional A(M) G (E/L)'^, that 

DjSif)\f. = A(M) -m, mifl,) = M, (8) 

or 

DfSif)\f,^ = ■ A(M), m{fl,) = M. (9) 
From equation Q we get, 

m(DM/M) = hE/L), (10) 

where we have indicated the space in which the unit operator is acting. Next, using the 
latter expression, we transform the differential of the equation (jSI), 

DMl^={m{D]S)-f}m^)-\ (11) 

and, consequently, 

DMfh = {DjS)j}y{m{DjS)j}yr\ (12) 

Notice that, elsewhere in equation (fT^ . operator {DjS)~^ acts on the hnear functionals 
from imm^. These functionals are precisely those which become zero on L (that is, on 
kerm), or, which is the same, those which can be represented as functionals of macroscopic 
variables. 

The tangent space to the quasiequilibrium manifold in the point flj is the image of 
the operator DmJm- 
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im {DMf*M) = iDjS)j}imm^ = {DjS)j}AnnL (13) 

where AnnL (the annulator of L) is the set of hnear functionals which become zero on L. 
Another way to write equation (fT3j) is the following: 

X e im {Dm flj) ^ {DjS)rJx,y) =0, yeL (14) 

This means that im (DmI'm) is the orthogonal completement of L in with respect to 
the scalar product, 

{x\y)f.^ = -iD}S)rJx,y). (15) 

The entropic scalar product (jl5|) appears often in the constructions below. (Usually, 
this becomes the scalar product indeed after the conservation laws are excluded). Let us 
denote as Tf^ = im^DMfli) the tangent space to the quasiequilibrium manifold in the 
point fl^. An important role in the construction of quasiequilibrium dynamics and its 
generalizations is played by the quasiequilibrium projector, an operator which projects E 
on Tf*^ parallel to L. This is the orthogonal projector with respect to the entropic scalar 
product, TTf*^ : E^Tf;y. 

TTf.^ = {DMrM)Mm= {D)S)-ly (m {Dp)-!^m^y'm. (16) 

It is straightforward to check the equality 7r|, = vr^.^, and the self-adjointness of nj*^ with 
respect to entropic scalar product (fT^. Thus, we have introduced the basic constructions: 
quasiequihbrium manifold, entropic scalar product, and quasiequilibrium projector (Fig. 

m 

2.2 Thermodynamic projector 

The construction of the quasiequilibrium allows for the following generalization: Almost 
every manifold can be represented as a set of minimizers of the entropy under linear 
constraints. However, in general, these linear constraints will depend on the point on the 
manifold. 

So, let the manifold Q = J'm C U he given. This is a parametric set of distribution 
function, however, now macroscopic variables M are not functionals on i? or f/ but just 
parameters defining the point on the manifold. The problem is how to extend definitions 
of M onto a neighborhood of J'm in such a way that Jm will appear as the solution to the 
variational problem: 

S{f) max, m(/) = M. (17) 

For each point /m, we identify Tm G E, the tangent space to the manifold Q in /m, 
and subspace Lm C E, which depends smoothly on M, and which has the property, 
Lm Tm = E. Let us define m(/) in the neighborhood of /m in such a way, that 
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QE Projector 




Kf = (D^S)"^m"^(m(D^S)"^m^)"^m 



M 

gure 3: Quasiequilibrium manifold Q, tangent space Tf*^Vt, quasiequilibrium projector 
• , and defect of invariance, A = Ap = J — TCf* (J). 
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m(/) = M, if / - fM e Lm. (18) 
The point Jm will be the solution of the quasiequilibrium problem ()17|) if and only if 



DfS{f)\f^^e Ann Lm. (19) 

That is, if and only if Lm C kerD/5'(/)|j^^. It is always possible to construct subspaces 
Lm with the properties just specified, at least locally, if the functional DfS{f)\j:^^ is not 
identically equal to zero on Tm- 

The construction just described allows to consider practically any manifold as a 
quasiequilibrium. This construction is required when one seeks the induced dynamics 
on a given manifold. Then the vector fields are projected on Tm parallel to Lm, and this 
preserves intact the basic properties of the quasiequilibrium approximations. 

It was proven [SEl 122] the theorem of uniqueness of the the thermodynamic 
projector: There exists the unique operator which transforms the arbitrary vector field 
equipped with the given Lyapunov function into a vector field with the same Lyapunov 
function (and also this happens on any manifold which is not tangent to the level of the 
Lyapunov function). 

Thermodynamic projector is constructed in the following way: Assume that the man- 
ifold Q G U is given, / G f2 and T is the tangent space to the manifold Q in the point 
/. Let us describe the construction of the thermodynamic projector onto tangent space 
T in the point /. 

Let us consider Tq that is a subspace of T and which is annulled by the differential S 
in the point /: 

To = {aeT\iDfS)ia) = 0} (20) 

If To = T, then the thermodynamic projector is the orthogonal projector on T with 
respect to the entropic scalar product (|)/ (fT^. Suppose that Tq ^ T. Let eg & T, eg -L Tq 
with respect to the entropic scalar product (|)/, and {DfS){eg) = 1. These conditions 
define vector eg uniquely. 

The projector onto T is defined by the formula 

PiJ) = PoiJ)+egiDfS)iJ) (21) 

where Pq is the orthogonal projector with respect to the entropic scalar product (|)/. For 
example, if T a finite-dimensional space, then the projector (j?T|) is constructed in the 
following way. Let basis in T, and for definiteness, {DfS){ei) ^ 0. 

1) Let us construct a system of vectors 

6i = Cj+i - AjCi, (i = 1, - 1), (22) 

where Aj = {DfS){ei+i)/{DfS){ei), and hence {DfS){bi) = 0. Thus, is a basis in 

To. 
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2) Let us orthogonalize {^i}"^^ with respect to the entropic scalar product (|)/ (fT3j). We 
have got an orthonormal with respect to (|)/ basis {gi}i~^ in Tq. 

3) We find eg & T from the conditions: 

(e, I g,)f = 0, (2 = 1, .., n - 1), (D/^)(e,) = 1. (23) 

and, finally we get 

n-l 
i=l 

If {DfS){T) = 0, then the projector P is simply the orthogonal projector with respect 
to the entropic scalar product. This is possible, for example, if / is the global maximum 
of entropy point (equilibrium). Then 

n 

i=l 

If {DfS){T) = and / is not equilibrium {Q is tangent to the to the level of the en- 
tropy), then the dynamic / = J{f) can be projected on Q with preservation of dissipation 
only if {DfS){J{f)) = in this point. 

2.3 Quasiequilibrium approximation 

Let a kinetic equation be defined in U: 

ft = ^^^^ 

(This can be the Liouville equation, the Boltzmann equation, and so on, dependent on 
which level of precision is taken for the microscopic description.) One seeks the dynamics 
of the macroscopic variables M. If we adopt the thesis that the solutions of the equation 
(EHjl of interest for us begin on the quasiequilibrium manifold, and stay close to it for 
all the later times, then, as the first approximation, we can take the quasiequilibrium 
approximation. It is constructed this way: We regard / as the quasiequilibrium, and 
write, 

^ = ^(^(/m))- (27) 

With this, the corresponding to M point on the quasiequilibrium manifold moves accord- 
ing to the following equation: 

df* 

= {DMfM)m{J{fl,)) = i^rjUli). (28) 

where ttj.^ is the quasiequilibrium projector (fTHj) . It is instructive to represent solutions 
to equations of the quasiequilibrium approximation in the following way: Let Tr{f) 
be the shift operator along the phase fiow of equation (j26p (that is, T^-lf) is solution to 
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equation ()26|) at the time t = t with the initial condition / at t = 0). Let us take the initial 
point /o = /mq, and set, /1/2 = T^(/o), Mi = m(/i/2), /i = /m^, . . ., /n+1/2 = T^(/n), 
M„+i = m{fn+i/2), ■ ■ ■■ The sequence /„ will be termed the Ehrenfest's chain. We set, 
frinr) = fn- Then, /r(t) fit), where fit) is the solution to the quasiequilibrium 
approximation (j^Hj) . as r ^ 0, n — 00, nr = t. 

Let us notice that the way the entropy evolves in time according to the Ehrenfests' 
chain is defined in the limit r ^ solely by the way it evolves along trajectories of the 
kinetic equation (pUj). Indeed, is the point of maximum of the entropy on the subspace 
defined by equation, m(/) = M. Therefore, for 

S{f n+1/2) - S{f n+i) = o(||/„+i/2 - /„+i||) = o(r), 

it holds 

J2\S(fn+l/2)-S{fn+l)\=o{nT)^0, 

n 

for T 0, n ^ 00, riT = const. This simple observation has a rather important im- 
plication: Let us denote as dSif)/dt the entropy production due to the original kinetic 
equation and as {dS{flf)/dt)i its derivative due to the quasiequilibrium system (j2HI)- 
Then, 

{dS{fl,)/dt), = dSif)/dt\^^^,^. (29) 

Let us give a different formulation of the latter identity. Let us term function S{M) = 
S{flf) the quasiequilibrium entropy. Let us denote as dSiM)/dt the derivative of the 
quasiequihbrium entropy due to the quasiequihbrium approximation (|77|). Then, 



dS{M) dSi^f) 



dt dt 



(30) 



From the identity ()29j] . it follows the theorem about preservation of the type of 
dynamics: 

(i) If, for the original kinetic equation (jSEI), dSif)/dt = at / = flj, then the entropy 
is conserved due to the quasiequilibrium system pHjl . 

(ii) If, for the original kinetic equation dSif)/dt > at / = flj, then, at the 
same points, /Xf, dS{M)/dt > due to the quasiequilibrium system (j2Zj). 

The theorem about the preservation of the type of dynamics^ demonstrates that, if 
there was no dissipation in the original system ()2fij) (if the entropy was conserved) then 
there is also no dissipation in the quasiequilibrium approximation. The passage to the 
quasiequilibrium does not introduce irreversibility (the reverse may happen, for example, 
there is no dissipation in the quasiequilibrium approximation for hydrodynamic variables 
as obtained from the Boltzmann kinetic equation; though dissipation is present in the 
Boltzmann equation, it occurs in different points but on the quasiequilibrium manifold of 

•^This is a rather old theorem, one of us had published this theorem in 1984 already as textbook 
material ([27j, chapter 3 "Quasiequilibrium and entropy maximum", p. 37, see also the paper j28j). but 
from time to time different particular cases of this theorem are continued to be published as new results. 



23 



Ehrenfestics 




M = F(M) 

Figure 4: Ehrenfest's chain over the quasiequihbrium manifold, and trajectory of the 
macroscopic dynamics, M = F{M). 
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local Maxwellians the entropy production is equal to zero). The same statement also hold 
for the thermodynamic projectors described in Section lT^ On the other hand, the entropy 
production in the quasiequilibrium state is the same, as for the quasiequilibrium system in 
the corresponding point, hence, if the initial system is dissipative, then quasiequilibrium 
entropy production is nonnegative. 

Usually, the original dynamics ()26|) does not leave the quasiequilibrium manifold in- 
variant, that is, vector field J{f) is not tangent to the quasiequilibrium manifold in all 
its points flf. In other words, the condition of invariance, 

{l-nf.JJ{f;,) = 0, (31) 

is not satisfied. The left hand side of the invariance condition ()31|) is of such an outstand- 
ing importance that it deserves a separate name. We call it the defect of invariance, and 
denote it as Aj*^. It is possible to consider the invariance condition as an equation, and 
to compute corrections to the quasiequilibrium approximation flj in such a way as to 
make it "more invariant" . In those cases where the original equation (j26|) is already dissi- 
pative, this route of corrections, supplemented by the construction of the thermodynamic 
projector as in Section UTR leads to an appropriate macroscopic kinetics ^Tj . 

However, here, we are mainly interested in the route "from the very beginning" , from 
conservative systems to dissipative. And here solving of the invariance equation does not 
help since it will lead us to, while "more invariant" , but still conservative dynamics. In all 
the approaches to this problem (passage from the conservative to the dissipative systems), 
dissipation is introduced in a more or less explicit fashion by various assumptions about 
the "short memory". The originating point of our constructions will be the absolutely 
transparent and explicit approach of Ehrenfests. 



3 Natural projector and models of nonequilibrium 
dynamics 

3.1 Natural projector 

So, let the original system be conservative, and thus, dS{f)/dt = 0. The idea of 
Ehrenfests is to supplement the dynamics (j26|) by "shakings". Shakings are external 
perturbations which are applied periodically with a fixed time interval r, and which lead 
to a "forgetting" of the small scale (nonequilibrium) details of the dynamics. For us 
here the shaking is the replacement of / with the quasiequilibrium distribution 
In the particular case which was originally considered in by Ehrenfests, the macroscopic 
variables m{f) were the averages of / over cells in the phase space, while fm{f) 
cell-homogeneous distribution with with the constant density within each cell equal to the 
corresponding cell- average of /. As we have already mentioned it, in the limit r — 0, one 
gets back the quasiequilibrium approximation - and the type of the dynamics is preserved. 
In this limit we obtain just the usual projection of the vector field J(/) on the tangent 
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bundle to the quasiequilibrium manifold. So, the natural question appears: What will 
happen, if we will not just send r to zero but will consider finite, and even large, r? In 
such an approach, not just the vector fields are projected but segments of trajectories. We 
shall term this way of projecting the natural. Let us now pose the problem of the natural 
projector formally. Let Tt{f) be the phase flow of the system (j^Bj) . We must derive a 
phase flow of the macroscopic system, Ot(M) (that is, the phase flow of the macroscopic 
system, dM/dt = F{M), which we are looking for), such that, for any M, 

miUfM)) = QriM). (32) 

That is, when moving along the macroscopic trajectory, after the time r we must obtain 
the same values of the macroscopic variables as if we were moving along the true micro- 
scopic trajectory for the same time r, starting with the quasiequilibrium initial condition 

(Fig. ED. 

It is instructive to remark that, at finite r, the entropy growth follows immediately 
from equation (j22I) because S{f) < S{f.^i^j^). The difference of the values of the entropy is 
of the order ||/ — /^(j) ||^, for the time r, thus, the first non- vanishing order in the entropy 
production will be of the order of r. Let us find it. 

We shall seek F in terms of a series in r. Let us expand F and both the sides of the 
equation (|32p in powers of r to second order, and find the expansion coefficients of F 

Trifo) = /o + df/dt\fj + d'f/de\^y/2) + o(r2), 
e.(Mo) = Mo + dM/dt\,,T + d''M/dt'\^y/2) + o{t^), 
df/dt\^^ = J(/o), rfVM'l/o = DfJ{f)\^^J{fo), 
dM/dt\^^ = F{Mo), d'M/dt%^ = DmF{M)\,,F{M,), 
F{M) = Fo{M) + tFi{M) + o(r). 

Using these expansions in the condition for natural projector ((221) , S^^i 

/o = /a/q) 

m(/o) +rm(J(/o)) + {r^ /2)D fJ{f)\^^J{U) + o{t^) 



Mo + rFo(Mo) + r^F,{M^) + {t^ /2)DmF{M)\j^,F{M^) + o(r 



2^ 



whereupon, 



Fo(M)=m(J(/l,)), 

Fi(M) = (1/2) {m{DfJ{f)\f.J{fl,)) - DmF,{M)\,,F,{M) 



"^In this well known work Lewis expanded only the right hand side of equation ^VZ\ . and did not 
do the same also with the left hand side. There were some hidden reason for this "inconsistency" : it 
was impossible to obtain the Boltzmann equation without such a deformation of expansion. We stress 
that our approach of matched expansion for exploring the coarse-graining condition is, in fact, the exact 
(formal) statement that the unknown macroscopic dynamics which causes the shift of M on the left 
hand side of equation (|32|l can be reconstructed order-by-order to any degree of accuracy, whereas the 
low-order truncations may be useful for certain physical situations. A thorough study of the cases beyond 
the lower-order truncations is of great importance which is left for future work. 
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Natural projector 



f = J(f) T,f(0) = f(t) 




©,f (0) = f (t) 



Figure 5: Projection of segments of trajectories: The microscopic motion above the man- 
ifold Q and the macroscopic motion on this manifold. If these motions began in the same 
point on Q, then, after time r, projection of the microscopic state onto Q should coin- 
cide with the result of the macroscopic motion on Q. For quasiequilibrium Q projector 
Tc : E -^n acts as Tc{f) = f^^^y 
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Thus, the approximation Fq is the quasiequihbrium, and using this fact in the expression 
for Fi, after some transformation, we derive, 

F, = (l/2){m(D;J(/)|^./(/;j)-DM(m(J(/;,)))m(J(/;,))} 

= (1/2) {m{DfJif)\^,Jifl,)) - miDfJif)\^.DMfM)miJifl,))} 

= (l/2)m {DfJif)\^,^ [Jin,) - DMfMm{J{fM))] 

= {l/2)m[DfJ{f)\^,Jl-irrJJ{rM: 

= {l/2)m{DfJU)\r,,^n 

Thus, the final form of the equation for the macroscopic variables M may be written: 

^ = F(M) = m(J(/;,)) + (r/2)m(D^J(/)|^. A^.J + o{t% (33) 



It is remarkable the appearance of the defect of invariance in the second term (propor- 
tional to r): If the quasiequihbrium manifold is invariant with respect to the microscopic 
dynamics, then F{M) is quasiequihbrium. 

Let us compute the production of the quasiequihbrium entropy S{M) = S{flj) due 
to macroscopic equations (IHHj). neglecting the higher-order term o(r^). 

dSir^)/dt = (r/2)D;^(/)|^,^vr^. D;J(/)|^. A^.^. 

We notice that, 

^/5(/)l,;//|,=/^/5(/)|,.^, 
because ttj*^ is a projector, and also because the thermodynamic condition 

ker TTj.^ C ker DjS(/)|^,^ 

which follows from the definition of quasiequihbrium (jHI). Next, by our assumption, the 
system (H^ conserves the entropy, 

dS{f)/dt = DfS{f)\fJ{f) = (}. 

Let us differentiate the latter identity: 

DjSUlljJif) + DfSUllfDfJUllf = 0. (34) 
Thus, due to the right hand side of equation 

= (r/2)D,5(/)|,,D,J(/)|,, A,.^ 



-(r/2) [Dp{f)\^,/{fl,)) A 



n 



M 
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where we have used notation for entropic scalar product ()15|) . Finally, 



A/I, = (1 - ^rJJUh) = (1 - ^rJJifu) 



whereas projector iXf*^ is self-adjoint in the entropic scalar product (fT^ . Thus, (J(/|j)| Aj.^)y.^ = 



Thus, the quasiequilibrium entropy increases due to equation of macroscopic dynamics 
(I33|) in those points of the quasiequilibrium manifold where the defect of invariance is not 
equal to zero. This way we see how the problem of the natural projector (projected are 
not vector fields but segments of trajectories) results in the dissipative equations. For 
specific examples see [22] where the second term in equation (j33j) results in viscous terms 
in the Navier-Stokes equations, diffusion and other dissipative contributions. However, 
it remains the undetermined coefficient r. Formula (j35p gives the entropy production 
just proportional to the time interval between subsequent coarse-graining. Of course, 
this could be true only for small enough r, whereas we are mostly interested in the limit 
T ^ oo. It is only in this limit where one can get rid of the arbitrariness in the choice of 
r present in equations (jH^ and In order to do this, we need to study more carefully 
the structure of the trajectories which begin on the quasiequilibrium manifold. 

3.2 One-dimensional model of nonequilibrium states 

In the background of many derivations of nonequilibrium kinetic equations there is present 
the following picture: Above each point of the quasiequilibrium manifold there is located 
a huge subspace of nonequilibrium distributions with the same values of the macroscopic 
variables, as in the quasiequilibrium point. It is as if the motion decomposes into two 
projections, above the point on the quasiequilibrium manifold, and in the projection on 
this manifold. The motion in each layer above the quasiequilibrium points is extremely 
complicated, but fast, and everything quickly settles in this fast motion. 

However, upon a more careful looking into the motions which begin in the quasiequi- 
librium points, we will observe that, above each point of the quasiequilibrium manifold 
it is located just a single curve, and all the nonequilibrium (not-quasiequilibrium) states 
which come into the game form just a one-dimensional manifold. It is namely this curve 
the construction of which we shall be dealing with in this section. 

This is the curve of the primitive macro scopically definable ensembles. These ensembles 
appear as the result (for t > 0) of motions which start from the quasiequilibrium state 
(at t = 0). 

For each value of the macroscopic variables M, and for each time r, we define M_t- by 
the following equality: 



(Ar |Ap )f* , and 




(35) 



(36) 
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In other words, M^r are those values of macroscopic variables which satisfy Bt-(M_t-) = M 
for the natural projector ((221) • Of course, it may well happen that such M.,- exists not 
for every pair (M, r) but we shall assume here that for every M there exists such tm > 
that there exists M_j. for < r < Tm- 

A set of distributions, qM,T = TT-{flj_^), forms precisely the desired curve of nonequi- 
librium states with given values of M. Notice that, for each r, it holds, m{qM,T) = M. 
The set {(Ja/.t} for all possible M and r is positive invariant: If the motion of the system 
starts on it at some time to, it stays on it also at t > to- If the dependence g^.r is known, 
equations of motion in the coordinate system (M, r) have a simple form: 

dr 
dM 

— = m{J{qM,r))- 

The simplest way to study qm^t is through a consideration of a sequence of its deriva- 
tives with respect to r at fixed M. The first derivative is readily written as. 



dqM,T 



dr 



r=0 



>/(/m)-^/x/(/m) = A/X,- (3^ 



By the construction of the quasiequilibrium manifold (we remind that L = ker m), for 
any x G L, 

Sifli + rx) = Sin,) - (rV2)(x|x);.^ + air'). 

Therefore, 

S{qM,r) = Sin,) - (rV2)(A,.JA,.;,.^ + o{t'). 

Thus, to first order in r, we have, as expected. 

qM,r = /m + rAj.^ +o(r). 

Let us find qM,T to the accuracy of the order o(r^). To this end, we expand all the 
functions in equation to the order of o(r^). With 

M_r = M- rmiJif*^)) + r'B{M) + o{t^), 
where function B is yet unknown, we write: 

/m_. = fli - tDmIIMJUIi)) + T'DMfMBiM) + (rV2)A2(M) + oir^), 

where 

df^ 



A,iM) = -^^^fj^^^^-" , (39) 

t=o 



and 
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Tt.{x + ra) = X + ra + tJ{x) + T'^D.j.J{x)\^a 
+ (rV2)Z},,J(x)LJ(x)+o(r2), 
TrifMj = rM-rDMflMJirM)) + r'DMfMB{M) + {T'/2)A2iM) 
+rJ{rM) - T'DjJ{f)\^,DMflMArM)) 
Hr'/2)DfJU)\rjUM) + o{r') 
= ft, + rAf., + (rV2)A2(M) + (rV2)D^ J(/)|^.^(l - 27r;.J J(/X,) 

+T^DMrMB{M)+0{T^). 

The latter somewhat lengthy expression simplifies significantly under the action of m. 
Indeed, 

m{A2{M)) = S[M + tm{J{rM))]/dt^ = 0, 
m{DMrM) = 1. 

Thus, 

HTrifliJ) = M+ ir'/2)miDfJif)\f,Jl - 27r;.J J(/;,)) + r'BiM) + (.(r^), 

BiM) = il/2)miDfJif)\f,j2nj,^^ - 1)J(/;,)). 
Accordingly, to second order in r, 

qM,r = T^ifhJ (40) 
= /l, + rA^;^ + (rV2)A2(M) 

+ (rV2)(l - 7r^.jD^J(/)|^./l - 27r^*J J(/;,) + oir'). 

Notice that, besides the dynamic contribution of the order of (the last term), there 
appears also the term A2 ((221) which is related to the curvature of the quasiequilibrium 
manifold along the quasiequilibrium trajectory. 

Let us address the behavior of the entropy production in the neighborhood of flj. 
Let X E L (that is, m{x) = 0). The production of the quasiequilibrium entropy, al.j{x), 
equals, by definition, 

aUx) = DmSU;,) ■ m{J{fl, + x)). (41) 

Equation ()4H) gives the rate of entropy change under the motion of the projection of the 
state onto the quasiequilibrium manifold if the true trajectory goes through the point 
flj + X. In order to compute the right hand side of equation (HH), we use essentially the 
same argument, as in the proof of the entropy production formula ()35|) . Namely, in the 
point /i^, we have L C kerD/S'(/)L* , and thus DfS{f)\f, ttj* = DfS{f)\f, . Using this, 

J ]\,T J A.'f A/ J A,I 
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and the fact that entropy production in the quasiequihbrium approximation is equal to 
zero, equation PT|) may be written, 

aUx) = DfSif)\ iJ{fl, + x) - JUm)). (42) 

To the hnear order in x, the latter expression reads: 

al,{x) = DfS{f)\^,pfJ{f)\j,^x. (43) 

Using the identity (jH^ . we obtain in equation (jlSj) . 

al,{x) = -Dp{f)\^,JJ{fl,),x) = (J(/l,)|x);.^. (44) 

Because x E L, we have (1 — tt/-* )x = x, and 

{■JifMMn, = (J(/M)|(l-7r;x.)^)/M 

Thus, finally, the entropy production in the formalism developed here, to the linear order 
reads, 

^lM = {^r,Mnr (45) 
3.3 Stability of quasiequilibrium manifolds 

The notion of stability does not cause essential difficulties when it goes about an invariant 
manifold, it is stable if, for any e > 0, there exist such 6 > that a motion which has 
started at t = at the distance (in some appropriate sense) less than 6 from the manifold 
will not go away further than e at any t > 0. 

However, this is not so for a non-invariant manifold, and, probably, it is not possible to 
give a useful for all the cases formalization of the notion of stability of the quasiequilibrium 
manifold, in the spirit of motions going not far away when started sufficiently close to 
the manifold (indeed, what is here "sufficiently close" and "not far"?). In spite of that, 
expression (jiUI) gives important opportunity to measure the stability. Indeed, let us 
consider how the entropy production depends on r, that is, let us study the function, 

aMir) = {Af,JqM,r)fi,. (46) 

It is natural to expect that aMi^) initially increases, and then it saturates to some limiting 
value. The question is, however, how function o"m(t) behaves at t = 0, is it concave or 
is it convex in this point? If function o"m(t) is concave, (i^crjvf (r)/(ir^|^^Q < 0, then the 
speed with which it grows reduces immediately, and one can even estimate the limiting 
value, 

cr^ = lim o-Af(r), 
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using the first Pade approximate: 

aM{r) = ar/(l + br) = ar - abr'^ + ... (47) 
* _ n_ 2(c;(TM(r)/^r|^^o)^ 

Concavity of o"Af(r) at r = {d'^o'M{T) /dr'^l^^^ < 0) is analogous to a soft instability: 
The motion does not run too far away, and it is possible to estimate where it will stop, 
see equation (jTfj) . However, if (i^o"M(T)/c?T^|^=o ^ 0, then this is analogous to a hard 
instability, and none of the estimates like ()47|) work. Thus, everything is defined by the 
sign of the scalar product. 



r=0 



If this expression is negative, then the Pade estimate (jlTj) gives: 

{^ri,\A,{M) + DjJU)\fAi-2i^n,)J{rM))n; ^ ^ 



In the opposite case, if the sign of the expression ()48|) is positive, we call the quasiequi- 
librium manifold unstable. 

Equation (jl^ allows us to estimate the parameter r in the equations of the method 
of natural projector. To this end, we make use of equation (jH^ : 

(^/2)(A/lJA/I.)/I.=</, 

whereupon, 

4(Ap lAp )f* 

_ _ , ^ J M ' J M I J M (50) 



(A/lJ^2(M) + DfJ{f)\^,^{l - 2nrJJirM))fl, ' 

if the denominator assumes negative values. In this case, there are no free parameters left 
in equation (j^Hjl . 

Above, the parameter r, or the time of "leaving the initial quasiequilibrium condition" , 
has been appearing explicitly in the equations. Except for the case of linear quasiequi- 
librium manifolds where the formal limit r ^ oo can be addressed to derive generalized 
fluctuation-dissipation relations [SI], this may be not the best way to do in the general, 
nonlinear case. 



3.4 Curvature and entropy production: Entropic circle and first 
kinetic equations 

In a consequent geometric approach to the problem of constructing the one-dimensional 
model of nonequilibrium states it is sufficient to consider the entropic parameter, SS = 
S*{M) — S. Within this parameterization of the one- dimensional curve of the nonequilib- 
rium states, one has to address functions o"m(A5'), rather than (lAf (r) (jlFIjl . whereas their 
Pade approximates can be constructed, in turn, from expansions in r. 
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In order to give an example here, we notice that the simplest geometric estimate 
amounts to approximating the trajectory qM,T with a second order curve. Given qM,T and 
(iM,T (001), we construct a tangent circle (in the entropic metrics, (|)/*^, since the entropy 
is the integral of motion of the original equations). For the radius of this circle we get, 

R = ^tdtiM^, (51) 

where 

qM,o = A/*^, 

■■ _ ■■ (gAf,o|A/»^)/»^Aj|^ 

q±M,o - qufl- 1^ V > 

qM,o = (1 - nrjDfJif)\f.Jl - 2nf,JJifl,) + {DM^rJ m{J{fl,)). 

Let us represent the microscopic motion as a circular motion along this entropic circle 
with constant velocity gM,o = A/*^. When the microscopic motion passed the quarter of 
the circle, the entropy production started to decrease and it became zero after the halve 
of the circle. Hence, after passing the quarter of the circle, this model should be changed. 
The time of the motion along the quarter of the model entropic circle is: 




\ J M' J M ' J M (52) 



After averaging on the 1/4 of this circle circular motion we obtain the macroscopic 



equations^ 



,n(j( n, + -fl^ + (l - -') R,^, ] I = MAfl,)) + (53) 



dt V V ^ l|A/;J| V ^/ M,o| 
2 R /„ \ f 2\ R 



71 \\Af* 



n 

where \\y\\ = sj {y\y) fi,- 

Equations ()53|) contain no undetermined parameters. This is the simplest example 
of the general macroscopic equations obtained by the natural projector. The coefficients 
(2/7r, etc.) can be corrected, but the form is more universal. The entropy production 
for equations is proportional both to the defect of invariance and to the radius of 
curvature: ^ 

aM = -R\\ArJ. (54) 

This equation demonstrates the thermodynamical sense of curvature of the curve of 
nonequilibrium states. The combination ci^yatur"^'^^^ dissipation. (It should 

be remained that all the scalar products and norms are entropic). 

""This averaging makes sense for conservative microdynamics, and for dissipative microdynamics the 
model of uniform circular motion along the entropic circle should be improved by taking into account the 
acceleration along the circle. 
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4 The film of non-equilibrium states 



4.1 Equations for the film 

The set qM,T in a space E generates a surface parameterized by "two variables" : A scalar, 
r > 0, and value of macroscopic variables, M, with condition 

M = m(gM,r). (55) 

We call this surface the film of non- equilibrium states or simply the film. It consists of the 
primitive macro scopically definable ensembles, the results (for t > 0) of motions which 
start from the quasiequilibrium state (at t = 0). 

For each r > the section of the film is defined: the set, qM,Ti for given r. It is 
parameterized by the value of M. For r = the section of the film coincides with the 
quasiequilibrium manifold. The film itself can be considered as a trajectory of motion 
of the section under variation of r G [0; +oo) (Fig. It is not difficult to write down 
equations of this motion using definition of qM,T '■ 

qM,T = Trf*M_^, (56) 

where Tr is shift in time in accordance with the original dynamical system, M_t- is defined 
with equation 
For small Ar 

qM,T+AT = qM-AM,T + J {qM,r)AT + o(Ar), (57) 
where AM = mJ{qM,T)A.T. Hence, 

= (1 - DMqM,Tm)J{qM,r)- (58) 
Initial condition for equation ()58|) is the quasiequilibrium 

gM,o = f*M- (59) 

Equation (jKH|l under initial condition defines the film of non-equilibrium states 
in the space E. This film is a minimal positive invariant set (i.e invariant with respect 
to shift Tj. by positive times r > 0), including the quasiequilibrium manifold, fl^. The 
kinetics we are interested in occurs only on this film. 

Investigation of non-equilibrium kinetics can be separated into two problems: 

i) Construction of the film of non-equilibrium states: solution of equation under 
initial conditions (f^. 

ii) Investigation of the motion of the system on the film. 

Naturally, it should be assumed that the film will be constructed only approximately. 
Therefore, the second problem should be separated in two again: 

iil) Construction of projection of initial vector field, J, on the approximately found film, 
and construction of equations for M and r. 
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Invariant Film 
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Figure 6: The film of nonequilibrium states as tlie trajectory of motion of the quasiequi- 
hbrium manifold due to microscopic dynamics. 
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ii2) Investigation and solution of equations for M and r. 

It should be especially emphasized that existence of the film is not significantly ques- 
tionable (though, of course, the deriving of the theorems of existence and uniqueness for 
()58|). can turn into a complicated mathematical problem). In contrast, existence 
of kinetic coefficients (viscosity etc.), and generally, the fast convergence of dM/dt to a 
certain dependence dM/dt{M) is essentially a hypothesis which is not always true. 

Below we will be solving, mainly, the problem of construction of equations: problems 
iil) and ii2). And we will begin with the problem ii2). Thus, let the film be approximately 
constructed. 



4.2 Thermodynamic projector on the film 



We need the projector in order to project the vector field on the tangent space. The 
idea of thermodynamic projector i20j consists of a description of every manifold (sub- 
ject to certain requirements of transversality) as the quasiequilibrium one. For this, one 
constructs a projection of a neighborhood of the manifold on it, and later, the required 
projector is obtained by linearization. 

The projection of the neighborhood on the manifold should satisfy essentially only 
one condition: a point of manifold must be the point of maximum of the entropy on its 
preimage. If the preimage of point /* is a domain in the affine subspace, Lf* G E, then 
required condition is: 



DfS*iLf,-n^O. 



(60) 



where Lf* — f* is already the linear subspace in E. 

For such projections, a dissipative vector field is projected into a dissipative one, and 
a conservative vector field (with the entropy conservation) is projected into a conservative 
one, i.e. the entropy balance is exact. Thus, let the film, qM,T, be defined. Let us construct 
for it a thermodynamic projector. 

Under small variation of variables M and r 



AqM,r = DMqM,rAM + DrqM,rAT + o{AM,AT), 

AS = DfS\g,,^^AqM,r + o{AM,AT). 

After simple transformations we obtain: 

1 



(61) 



At 



^f^\qM,T^rqM,T 

DTqM,TDfS\q^j 



o(AM, A5), 

DMqM,r/^M 



+ 



^ f^\qM,T^TqM,T _ 

^ DrqM,r^S + o{AM,AS). 



DfS\gj^j^D^qM,r 

From this formulae we obtain thermodynamic projector for J, TTta: 



J 



DrqM,TE>fS\q^j^ 

DfS\g^j^^DrqM,T 



DrqM,TDfS\g^,j^ 
DMqM,TmJ + —— -J. 

^f^\qM,r^TqM,T 



(62) 



(63) 
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For conservative systems the second term in ()63|) vanishes and we obtain: 

DT(lM,TDfS\gj^j. 



TTtd 



J 



DfS\y,j^^DrqM. 
The equation for M corresponding to (j64j) has the form: 

DrqAI,TDfS\qj^j 



DMqM,rmJ. 



(64) 



m 



DfS\q„l^DrqM,' 

It should be supplemented with the equation for 5*: 

dS 



DMqM,TmJ = mJ {qM,r] 



(65) 



dt 



0, 



(66) 



or for r, in accordance with (j^ . 

dr S - DfS\q^j^DMqM,rM D jSlq^^^Duqu^rM 



dt 



^f^\qM,T^TqM,T 



^f^\qM,T^TqM,T 



(67) 



where M is defined in accordance with (jUH)). The numerator in (jU7j) has a simple meaning: 
it is the rate of the entropy production by dynamic equations fl65p when r is constant (for 
frozen r). Expression ()67|) can be obtained from the condition of the constant entropy for 
the motion on the film in accordance with ()65I67|) . Equations ()65I67|) describe dynamics 
on the film (Fig. Ej). 

Let us further assume that condition (IHHll is satisfied: 



Qm,- 



/m + rAf* + o(r) 



In expressions ()62|65|67|) the denominator, D fS\q^j ^DrqM^r^ is present. For r ^ this 
expression vanishes: 



-Dr'?M,T|r=0 

DfS\f=rx 



Af* , 
0, 



for 



X G ker m, 
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m{Af*^) = 0, therefore D fS\q^j ^D^qM^r — > for r — > 0. For r ^ indeterminate forms 
0/0 appear in expressions (j62ll64l66l67|) . Let us resolve the indeterminate forms and 
calculate the corresponding limits. 

Two indeterminate forms are present: 



{DrqM,r){DfS\qj,jjDMqM,TmJ 

^fS\qM,T^TqM,T 



(69) 



and right hand side of equation (jU7|) . Let us evaluate the indeterminate form ()69p with 
the L'Hopital rule. We obtain: 



N^(T) 



^nDjS\r^^7rr^DjJif)\,. 



(Ap \ Af* )f* 
38 



(70) 



Dynamic Equation on The Film 




using identity similar to we obtain: 

Ap (Ap \Af* )f* 
- (A,. |A,. = 

In such a way, for r ^ 

TTtdl.M,. J(gM,.) - DMrumJifl,) + A^;^ = TT^;, J(/;,) + (1 - l^r^JUu) = J Ull) ■ (71) 

Similarly, after simple calculations we obtain that: 

dr ^ , , 

- ^ 1. for . ^ 0, (72) 

The fact that for r — the action of the thermodynamic projector on J becomes 
trivial, n^dJ = J, can be obtained (without calculations) from the construction of gM,T 
in vicinity of zero. We have chosen this dependence in such a way that J{qM,T) becomes 
transverse to the film for r — 0. This follows from the condition ()38|) . Let us emphasise, 
however, that derivation of the formulas ()64II67|) themselves was not based on ()38p. and 
they are applicable to any ansatz, gA/,r, not necessarily with the right behavior near the 
quasiequilibrium (if one needs such ansatzes for anything). 



4.3 Fixed points and "right asymptotics" for the film equation 

What is the dynamics of the film in accordance with equation (j3Hj) ? A naive expectation 
that qM,T tends to the stable point of equation leads to strange consequences. Stable 
point is the invariant manifold Qm- On this manifold 

J(gA/) = DuqumJiqu), (73) 

i.e. the projection of the vector field, J, onto qu coincides with J. Were the condition 
qM,T qm satisfied for r — >^ cx), the dynamics would become more and more conservative. 
On the limit manifold qu, the entropy should be conserved. This could lead to unusual 
consequences. The first of them is limited extendability backwards "in the entropy" . 

Let us consider the set of points for given M. Because of the existence of the 
limit, Tt-M_t — i> qM, for r — 0, the difference, S{M) — S{M^^) = ASr, is bounded on the 
half-axis, r e [0; +oo) : ASr < AS'oo(M). this means that it is impossible to get into the 
values of macroscopic variables, M, from the quasiequilibrium initial conditions. Mi, for 
that S{M) — S{Mi) > AS'oo(M). Assuming additionally a smoothness of qM and M_t-, we 
see that it is impossible to get into e— neighborhood of the quasiequilibrium state, M*, 
(over macro- variables) from the outside, from the quasiequilibrium initial conditions Mq, 
if S{Mq) < Se, where is e dependent threshold of the entropy. Thus, possible stable 
points of the equation (jSHI), regardless of their obvious interest, likely demonstrate exotic 
possibilities. The following "right asymptotics" correspond to our quahtative expectations 
for large r. Namely, it is expected that for the quite large r, M becomes, within a good 
precision, a function of M, and later does not depend on r : 

m{J{qM,r)) ^ M(M); (74) 
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with the entropy production: 

a{qM,r) = DMS{M)mJ{qM,r) -> a(M) > 0, (75) 

and, correspondingly, S{qM,T) —oo, t ^ oo. 

Already simple examples (linear in J) demonstrate that it is not so simple to construct 
such an asymptotic. Moreover, for reasonably built systems it probably does not exist. 
Indeed, let J(g) = Lq, we search for "the right asymptotic" in the form qM,T = ct{M) + 
Tb{M) +o(l). We obtain: 

mb = mLb = 0, 
ma{M) = M, 
Lb{M) - DMb{M)mL{a{M)) = 0, 

La{M) - DMa{M)mL{a{M)) = b{M). (76) 

Acting with the operator niL on the first equation, we obtain mL'^b{M) = 0; further, 
acting with the operator mL^, we obtain mL%{M) = 0, and so on. 
Thus, 

oo 

b{M) C Pi ker mL'^ = H. (77) 

fc=0 

Space H is L— invariant, therefore, it is possible to pass from the initial dynamics, / = Lf, 
to the dynamics in the factor-space. This does not change the dynamics of macroscopic 
variables because of the definition of H (f77j) . 

In such a way, instead of the right asymptotic equations, (fTUj) leads us again to the 
equation of the invariant manifold (6 = 0, a{M) determines the invariant manifolds.) 

4.4 Coarse-graining projector 

A construction of an exact projection of the microscopic dynamics on the macroscopic 
is meaningless, it has meaning only as an intermediate result. Really, generically, such a 
projection (finite segment of the trajectory M{t)) contains practically all the information 
about the Liouville equation. This is a bit too much. 

Moreover, there are no invariant manifolds with the dissipative dynamics for the finite- 
dimensional conservative systems. The conclusion is: every time explicitly, or sometimes 
implicitly, coarse-graining, or replacing of the system with something different, takes 
place. 

For example, there is no invariant manifold for the Liouville equation parameterized 
with the one-particle distribution function with dissipative dynamics on this manifold. 
The derivation of the Boltzmann equation requires some limit transitions. 

A few ways of coarse-graining are known, but essentially only two exist. The first one 
is related to distinguishing a manifold, M, and a projector, 11, on it; the manifold, M, 
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with the projector, 11, separate the "microscopic" (kerll) from the macroscopic, M. In 
particular, m(/) = m(n/). Then, a "new microscopic dynamic" is parameterized: instead 
of M = m{J{f)) the equahty 

M = m(J(n/)) (78) 

is used. 

Equation ()78|) determines, for example, quite different system of relations in the chain 
of derivatives M, M, M (original system can contain no such relations). For construction 
of "right asymptotics" relation ()76p is much more appropriate then the correct initial 
relation: for example, instead of mL'' in ()77p one can get m(Ln)^, and the kernel (LII)^ 
always contains the kernel 11. 

If the coarse-graining projector, 11, and manifold, M, are postulated, then the hypoth- 
esis of existence of the thermodynamic limit consists of existence of the limit: 

nqM,r-^fti forr^oo. (79) 

Nevertheless, one should not expect a precise fulfillment of equality (f7^ . this is still 
the case of exact projection (on M). For t ^ oo one can expect realistically only the 
smallness of the remainder: 

6{T,M,N) = UqM,r-ft, (80) 
and its elimination for N oo, where is number of particles. 

For a more precise formalization of this condition one should estimate S under r ^ oo, 
for example, such as: 

e{M, N) = l\\6{T,M,NWdT (81) 

and, further, investigate the thermodynamic limit N ^ oo. 

So far, however, the number of particles was not taken into consideration, and we dealt 
with only one fixed system. In this case we can't help but to assume that the value of 
b{M) (j81|) is "sufficiently small" . We notice that the problem of determining the dynamic 
for the thermodynamic limit of infinite systems is a very difficult problem. For infinite 
system even determining the energy, entropy, and other characteristics is not clear. We 
are talking not about the limit does not exist (the number of particles is always finite), 
but about the asymptotics for large N, therefore, strictly speaking, one needs not only to 
know the limits, but to estimate the reminding terms too. 

The second coarse-graining method consists of a decomposition of the system into 
small subsystems, and introduction of two incomparable time scales: micro and macro. 
The main assumption is that during an arbitrarily small macroscopic time period a small 
part of the system passes the micro-evolution within an infinitely long time. This leads to 
a quasi-chemical description: within each period a number of elementary atomic processes 
(events) takes place. For example, the derivation of the Boltzmann equation in frames of 
all formalisms is, in fact, reduced to this. We will return to considering this approach, 
but for now we'll concern ourselves with the coarse-graining projector for the film. 
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4.5 Choice of the coarse-graining projector, and layer-by- layer 
linearization 

The simplest choice of the coarse-graining projector is 

n(/) = (M(/), M(/)) = (m/, mJ(/)). (82) 

For many problems, for example, to investigate the invariance defect fl78p it is not 
necessary to place the manifold M in the initial space it is sufficient to investigate 
Mt. 

In the cases when one needs, after all, to have corresponding elements of E, a good 
choice could be the quasiequilibrium manifold corresponding to 11. The quasiequilibrium 
manifolds have an evident but important property. Let 

/ Ml ^ M2, / M2, (83) 

be a sequence of linear mappings, where mi, m2 are mappings "on", their images are 
whole corresponding spaces. 

Let, furthermore. Mi C -E be a quasiequilibrium manifold in E, corresponding to mi, 
M2 be a quasiequilibrium manifold of macro- variables, M2, corresponding to m2, M21 be 
a quasiequilibrium manifold in E corresponding to m2mi. Then 

mi(M2i)=m2, or M21 = m^i(M2i). (84) 

For the transition to the quasiequilibrium approximation this property reads simpler: 

U2U1 = f/21, (85) 

where Ui the corresponding to m^ procedure of the taking of the quasiequilibrium approx- 
imation. 

For each M2 both the point of the quasiequilibrium, Mj*(M2), and the linear manifold, 
m^^(M2), containing this point are defined. For each Mi the quasiequilibrium, /j'(Mi), 
and the linear manifold, m^^{Mi) containing this point are defined. As well (M2), and 
containing them (m2mi)~^(M2) are defined. 

Relations: 

/2*(M2) = /r(M;(M2)), 
m^\m^\M2)) = (m2mi)-i(M2). (86) 

are fulfilled. 

The quasiequilibrium manifold, /|(M2) C E parameterized by M2 lies on the quasiequi- 
librium manifold, f*{Mi) C E parameterized by Mi. For each M2 the set 

{/*(Mi)|/*(Mi) C (m2mi)-i(M2)} (87) 

forms the quasiequilibrium manifold in (m2mi)~^(M2) with the set of macroscopic vari- 
ables, m2 "^(M2), and the same entropy. For the projector, n(/) = (M(/), M(/)), it means 
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that for each M in the hnear manifold, on that m{f) = M, the quasiequihbrium manifold 
corresponding to the macroscopic variables M{f) = mJ{f) (if J(/) is a linear mapping 
on this manifold) is defined. 

The last remark leads us to an important construction named by us "layer-by-layer 
linearization". The filed J{f) could be presented in the form: 

The "layer-by-layer quadratic entropy" has special importance for the theory of non- 
linear equations (|SS|) : 

SlU) = - (i/2)(/ - - r^if))r^^^y (89) 

Let us remind that the bilinear form, is generated by the negative second differ- 

ential of the entropy at the point /. 

The layer-by-layer linearized equations allowed us to add more moment equations and 
construct the quasiequihbrium approximations using the entropy (jH^ . This is especially 
important for the moments which are time derivatives M, M and so on. 

Application of the layer-by- layer linearized equations (jH^ together with the layer- 
by-layer quadratic entropy ()89|) allowed us to construct a thermodynamically consistent 
theory of the moment equations for the Boltzmann equation [32 , 33 . 

It is convenient to supplement the quasiequihbrium, fl^, with the quasiequihbrium for 
additional macro-variables M, 

Mif) = miDfJif)\f.^^J), (90) 

in two stages: i) supplementing by the entropy production, ii) and later by the conserving 
part of the entropy. 

i) We supplement M by the entropy production. In the layer-by-layer linear approxi- 
mation 

^(/) = (^/*(,)l/-A/))/,;(,) (91) 

(as was already determined, see ()45|) ). The quasiequihbrium manifold corresponding to a 
in the layer over fl^ has the form: 

/m.. = /m+ • (92) 

(^/mI^/m)/? 



I M 



Quasiequihbrium projector in the layer is: 



(At* \Af* ) 

\ Jm' JM' f 



M 



ii) We distinguish in L = DJ{f)\f*^ the conservative (conserving the entropy) part 
over fli : 
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This corresponds to the situation when we have fixed = {Af*^\(f), and we consider 
the motion in the layer for fixed a. For this: 



Lm^ - L^f^ + ^ ^ ^ - L,,^ + ^ I ^ . (93) 

The quasiequihbrium manifold corresponding to L^j in the layer over flj could be 
constructed in the following way: we search for the kernel of L^j, in L (the set of all 
solutions to equations L^^p = 0, rmp = 0). We define it as K. The orthogonal complement, 
K-^, to K in the scalar product, (|)/*^, is the corresponding manifold. For each point from 
the image, L^, on L, G L^^(L), there exists unique ip G Kf, such that LmV = V'- We 
define it as = {L'^^)~^ {ip) . As a result, for every ip G L1j{L) 

fk<^,^ = fli + ,^ . ^ + (^m)"V. (94) 

The second and third terms in ()94|) are reciprocally orthogonal in the scalar product 

(l)/.v 

4.6 The failure of the simplest Galerkin approximations for con- 
servative systems 

The simplest approach to the problem is connected to the Galerkin approaches: one con- 
siders a projection of the vector field, J{f), onto the manifold in question and investigates 
the obtained motion equations. It is not difficult to make sure that for conservative sys- 
tems such an approach is unfruitful. If the orthogonal projection, (|)/*^, is taken, then 
in the linear within the layer approximation only quasiequihbrium approximations with 
increased number of moments could be obtained. For the dissipative systems, in con- 
trast, such a way leads to quite satisfactory results. Thus, if for the Boltzmann equation 
and the hydrodynamic moments the invariant manifold is to be searched in the form 
fti ~ fti + '^(^)^/m' then we obtain the Navier-Stokes equations with the viscosity cal- 
culated within the first Sonine polynomials approximation. Using another scalar product 
simply leads to unphysical results. 

In order to specify appearing problems, let us give an example with a linear field, 
J(/) = Af, and quadratic entropy, S{f) = (l/2)(/|/). The conservativity of J means 
that for each / 

{f\Af) = (95) 

is fulfilled. 

The quasiequihbrium subspace corresponding to the moments M = mf is the orthog- 
onal complement, kerM. The quasiequihbrium projector, vr, is an orthogonal projector 
on this subspace, and does not depend on the point. For the defect of invariance Aj^ we 
obtain: 

A;. ={A- 7rA)fl,. (96) 
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Under Galerkin approximation we write 



qM,r = rM + a{M,T)Af*^. (97) 
Projector of the vector field on Af*^ is 



\ JAfI JM' 



(9^ 



Thus, we pass from the equation of motion of the film ()58|) to the Galerkin approxi- 
mation for a{M, t). 

d = l+ (99) 

— a — j— r'^ [UMajiTi- 





lAnAAf* ) 




\Af* ) 



(Ap lAp ) (Ap \Af* ) (Ap lAp ) ' ' (Ap lAp ) 

\ JM' JM' \ JM' JM' \ Jm' JM' \ Jm' JM' 



One can try to find the fixed points (solving a = 0). This is the projected invariance 
equation. Due to the properties of the operator A, and the self-adjoint projector, vr, we 
obtain for conservative systems 

{ArJAArJ = 0, (100) 



On the other hand, for the dissipative systems the form (jlOOp is negatively definite, 
and it is this form that determines the Navier-Stokes equations (in the first Sonine's poly- 
nomials approximation) for derivation of these equations from the Boltzmann equation. 
For the conservative equations this main part vanishes, and the second term in equation 
fllOOp . generally speaking, is sign- indefinite. 

The failure of the Galerkin approximations is even more obvious in the equations of 
motions on the film. Here everything is very simple: 

(Ap MAp ) 

\ J M ' J M' 

For the dissipative systems under frozen M, a relaxes to the stable point 



(Ap UAp )■ 

\ jm' Jm' 



(103) 



This fixed point is "the leading order term" in the solution of the invariance equation, 
d = (ITnn|) . 

For the conservative systems d = 1. This result was evident beforehand from the 
entropy production formula (jH^j) . and 

- S{f) = (l/2)(/|/) = (l/2)(7r/|vr/) + (1/2)((1 - 7r)/|(l - n)f). 
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4.7 Possible ways beyond the simplest Galerkin approximations 

The first way is an application of tlie projection operators metfiods Tlie film equation 
is considered for two sets of "variables": slow "macro- variables" , Ilg, and rapid 
"micro- variables" , (1 — Il)q (where 11 is a coarse-graining projector, see subsection I4.4|) . 

Next the rapid variables are eliminated, and the equation with retardation for the slow 
variables is written. This formally exact equation becomes tractable only after a sequence 
of additional approximations ("short memory", "Markovian models" etc.). The method 
is applicable to linear (linear within the layer) vector fields, J{q) = LmQ- The main 
problem is the computation of the coefficients including averaging along the trajectories 
of the rapid motion. 

The second way is an introduction of the dissipative part (using the thermodynamic 
limit) into the vector field, J{qM,T)- One adds into either a "relaxation" operator 



or operator 7P, simulating a random process. For example, if / is a function on a space 
X, then the typical form of P with the "detailed equilibrium" is fl^ 



with non-negative kernel Q{x, x') > 0, J Q{x, x')dx' = 0. 

As a result, the system becomes dissipative, and one can construct for it invariant 
manifolds that are stationary solutions for the film equation. They could be found either as 
a sequence ESI o'^; more effectively, based on the Newton method with the incomplete 
linearization [711 12H] • It becomes possible to use the Galerkin approximations, and so on. 

After this one makes the transition to the thermodjTiamic limit. It is suggested that 
the thermodynamic limit exists for the found invariant manifold, qffi'y), and, if later 7 
tends to zero, that a finite limit, Ilqf^{'y) exists. This limit is suggested for the definition 
of the macroscopic variables M. 

In some problems of dissipative kinetics (namely, in the problem of initial layer for 
the Boltzmann equation) it was found to be effective to approximate the trajectories by 
segments (with further smoothing and corrections, or without them). These segments 
were constructed in the following way: the initial direction of motion was taken, and / 
evolved along this direction for as long as it was possible to conserve the smooth evolution 
of the entropy. Further, the procedure was repeated from the obtained point (for details 



Unfortunately, in the problem of the initial layer for the conservative systems there are 
no stop points during the motion along the straight line (more precisely, the start of the 
motion itself can be considered as a stop point because under the linear approximation 
the relation (jlUUj) is valid). In the initial layer for the dissipative systems the motion of 
the system along the straight line x = tA in any case increases the entropy. For the 
conservative systems one needs to "rotate the phase" , and the models of motion are to be 



(104) 




(105) 



see uniiini)- 
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arcs of ellipses (in linear space), or the constant entropy lines, instead of the straight lines. 
In the film problem, as even the simplest examples show, the simplest good model is a 
general conic section. A simple example: J{f) = Af, A is generator of rotation around 
the axis with direction f = + aCy, M = x, the film is the lateral surface of the cone, 
obtained by rotation of the quasiequilibrium manifold, the axis {xe^}, around the axis 
{(pr}. For a < 1 the curve qM,T is an ellipse, for a > 1 it is a hyperbole, for a = 1 it is a 
parable. 



4.8 The film: Second order Kepler models 

The curve qM,T is a section of two manifolds: one of them is the result of motion of 
the quasiequilibrium manifold along the vector field J{f), and another one is the linear 
manifold flj + ker m. 

Already in the finite- dimensional space, and under linear approximation (J is linear, S 
is quadratic) we have an interesting geometrical picture: quasiequilibrium manifold is an 
orthogonal complement to kerm, A is the rotation generator, (kerm)-*- is rotated under 
action of e^"^, unknown curve is the section: 

(Z;^ + kerm) Pi e^^+ (kerm) ^, (106) 

where i?+ = [0; oo), flj e (kerm)-*-. 

The simplest model motion is a second order curve. However, it is not sufficient to 
know the first and the second derivatives. We need information about the third derivative. 
If we consider the curve qM,T as a trajectory in the Kepler problem, then the location, r, 
of the center of attraction (repulsion) is (Fig. |SI): 

r = qo- q , . , , 107 

where tq is the initial point where all the derivatives are taken. The force is: 

r - g 

r = 01 ; ■ 

(r — q\r — q)'^/"^ ' 

^ = (a\a? ... 

a > (attraction) if {q < 0; nna^ 
a < (repulsion) if (g > 0. 

It is necessary to specifically note that the Kepler problem defines an approximation 
of the trajectory ^m.t, but not the dependence on r. 

An important question is the finiteness of the film. Is the modeling motion finite? 
The answer is simple in terms of the Kepler problem 

||(7p a 



a' = {q\q){r-q\r-qf = {q\qf^-^^]^; (108) 



\r - qo\ 
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Second order models 



1 . The circle 



<q|q> 



Kepler models 



(qi I qi> 




r 



Figure 8: The definition of tlie second-order models. 
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or 



< 1. 



Here 



2||g±P||gP 

is the norm in the entropic scalar product, as it is usual. 



(110) 



4.9 Minimal second order models: entropic parable, and en- 
tropic circle 

In accordance with the film equation (j^Hj) . the following derivatives 



Qm,t 

(iM,T 



dqM,r/dT; 
d^qM,r/dT'; 



contribute to the construction of the second order Kepler models. 

There is a rougher construction leading to the two distinguished simplest second order 
models that uses only two derivatives. One of them is finite (the entropic circle), another 
one is infinite (the entropic parable). Both could be constructed for every point of the 
film (if q± ^ 0, otherwise, all second order models turn into straight lines). The circle 
was already used by us in subsection 13.31 in order to estimate stabilizating value of a. Let 
us remind that: 



R 



i'il'i) 



\/{q±\q±) 



[111) 



where q = dqM,T/dT, q = d'^qM^r/dr'^ 



(lL = q 



q{q\q)_ 



(I) is the entropic scalar product corresponding to the extension of the entropy at the 
point qM,To, or, for the linear (as well as linear within the layer) systems, at the point /|^. 
The concentric motion could be presented as: 



Qm,- 




cos 



M 
R 



T 



^0, 



The parable could be constructed simpler: 

gM,r - gM,ro = " ^o) + (l/2)g±(T " Tq)'^ , 



[112) 
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or even as: 



qM,r - qM,ro = ^(^ " ^o) + (l/2)g(r - Tof, (114) 

The difference between ()113|) and pi4j) is this: in the formula ()113|) the angle between 
d'iM,Tld'T\e=To and dgM^r/dr tends to n/2 for r — > oo, in the model ()114|) this occurs too 
if (gig) ^ 0. 

Note. In expressions ()112II114|) r coincides with the true r of the motion on the film 
only in the zeroth and the first orders. For the further calculations it could be necessary 
to recalculate r using true values of g in the projection on the trajectory. This is discussed 
below. 



4.10 The finite models: termination at the horizon points 

In order to construct a step-by-step approximation it is necessary to be able to solve two 
problems: the choice of the direction of the next step, and the choice of the value of this 
step. 

If the motion qM,T is along the straight line (dissipative systems), the direction of the 
next step is Qm^to (let us remind that gA/,ro is the defect of the invariance of the manifold 
Qm = (lM,To under fixed r = tq), and the value of the step should be taken in the direction 
to the stable point: to the point where direction of qM,T becomes orthogonal to initial 
one, qM,To (Fig- El)- Naturally, the current direction of qM,T is calculated with but 
approximately, with the frozen projector {DMqM,To instead of DmqM,Tf^)- 

For the conservative systems we have chosen the second order models instead of the 
linear ones. For finiteness of the models we need to define the moments of stop. It is 
suggested to operate in a manner similar to the case of the dissipative systems: to stop 
at the moment when the direction of motion is orthogonal to the initial one. In this case 
we will take the direction of motion along the model. 

Thus, if qM,To is a starting point of motion, and qM,To+e is a motion on the finite second 
order model, then condition for the transition to the next model is 

(^M,J%±^) = (115) 

(in the entropic scalar product). 

Let us call the horizon points such points, gM,ro+6»o5 where the scalar product (jllSp for 
the first time becomes zero (for < < this scalar product is positive). This notion is 
motivated by the fact that for 9 > 9q the motion on the second order model "disappears 
behind the horizon", and its orthogonal projection on the line parallel to qM,TQ starts to 
move back passing the same points for the second time. 

The convention about the change of the model in the horizon points seemed quite 
natural. The following sequence of calculations suggests itself (Fig. ITUI) : 

1) we pose that gM,o = /m; 

2) we calculate qufl, QMfi, ... in accordance with equation (50); 
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The Film for dissipative systems: 

step by step 



Terminal points 
in directions 




M2 M3 M4 M 



But ... there are no terminal points 
in directions A for conservative systems 

Figure 9: The stepwise construction of the film for dissipative system. First-order models: 
The motion along the defect of invar iance. 
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3) we construct the (finite) second order models, qufi] 

4) we find the horizon points, qM,9o{M)', 

5) then we take the manifold of the horizon points as a new initial manifold, and so 

on. 

At the first glance, this sequence contradicts the original statement of the film problem. 
The manifold qM,eo(M) does not have the form of qM,T for a fixed r and is not a shift of 
the quasiequilibrium manifold by the given time along the true microscopic equations of 
motion. 

The second difficulty was already mentioned: the time of motion along the modeling 
curve does not coincide with the true time, r. More precisely, it coincides only within the 
second order. However, global, not local approximation are consructed. Therefore, global 
corrections to the time, or ways to circumvent these corrections, are required. 

The following two sections are devoted to the elimination of these difficulties. 



4.11 The transversal restart lemma 

Let qM,T ij € [0; +oo)) be the solution to (jHHl) under initial condition (fS^ (the film). We 
call the transverse section of the film, qM,T, the manifold, qM,e{M), where 9{M) is a smooth 
function < e{M) <t<oo. 

Let the transverse condition be fulfilled. Namely, for every finite patch, M, that does 
not exclude equilibrium exists, e > such that in this patch 

II J iqM,e{M)) — DMqM,e{M)''TT'J iqM,e{M))\\ , . 

\\J{qMMM))\\ ^ ' 

in an appropriate (for example, entropic) norm. Let qM,T be the solution to under 
the initial condition qMfl = qivi,e{M)- Then the following transverse restart lemma is valid: 

(117) 

here qM\a;b] = {gA/,rk e [«;^]}- 

In order to prove^ this lemma, we notice that it is equivalent to the following statement. 
For every M the segment of the trajectory, Trfli (r G [0; t]), crosses the manifold qM,e{M), 
and only once. 

In order to demonstrate the unicity of the section, we consider the film in another 
coordinates, for each point q we set M and t,: q = T^/X/- 

In these coordinates the transverse condition excludes pleats on qM,e{M)- 
In order to demonstrate the existence of the segment T^/X/ (t G [0;t]) for the cross 
point, q*, with qM,e{M), we define in the neighborhood of the point flj on the quasiequi- 
librium manifold the mapping into the neighborhood of this section point. Image of the 
point flj is section of the trajectory T^flf (r G [0;t]) with the manifold qM,e{M) in the 

^Let us remind that in the degree of generality being used there are no proofs to the theorems of 
existence and uniqueness 
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The Film for conservative systems: 
second order steps 




M M 



Figure 10: The stepwise construction of the film for dissipative system. Finite second- 
order models: The motion starts in the direction of the defect of invariance, and stops 
when the direction of motion becomes orthogonal to the defect of invariance. 
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neighborhood of q* . Due to the transverse condition, it performs an isomorphism of the 
neighborhoods. Therefore, the set of M for which the section of the trajectory with 
<iM,e{M) exists is open. Furthermore, it is closed, because the hmit of section points is a 
section point (and segment [0;t] is compact). Obviously, it is not empty. Consequently, 
it is the set of all possible M. 

4.12 The time replacement, and the invariance of the thermo- 
dynamic projector 

Let the film be constructed as qufi-, where relation between 9 and r is unknown; r = 
t(M, 9), 9 = 9{M, 9), in order to determine this functions one needs to solve the equation 
obtaining from (jH^ with substitution g^.r = qM,e{M,T) (and projection, therefore, q is only 
an approximation). The calculation itself does not contain principal difficulties. However, 
a question arises: is it possible to escape the reverse replacing of time for the derivation of 
the kinetic equations? Another words, could we use the constructed geometrical object, 
the film, without an exact reconstruction of the time, r, on it? 

For positive answer to this question it is sufficient to demonstrate that the equations 
of motion, constructed with the thermodynamic projector (|65II67|) . describes the same 
motion on the film after the time replacement. 

This property of the thermodynamic projector is evident: deriving equations ()65H67|) . 
we did not use that r is the "true time" from the equation (j^ . and made the local 
replacement of variables, passing from AM, Ar to AM, A^*. 

In such a way, the thermodynamic projector is invariant with respect to the time 
replacement, and, constructing equations of motion, it is not necessary to restore the 
"true time". 

Results of this, and previous sections allow to apply the sequence of operations sug- 
gested in subsection 14.101 

4.13 Correction to the infinite models 

Let an infinite model qM,e, {9 G [0; +oo)), qMfi = be constructed for the film. Actually, 
it means that an approximation is constructed for the whole film qM,T (but not for its 
initial fragment, as it was for the finite models). Naturally, there arises a problem of 
correction to this initial approximation, and, in general, construction of a step-by-step 
computational procedure. 

The thermodynamic projector on the film is defined (jGlj) . Correspondingly, the invari- 
ance defect of the film is determined too 



AgAf.e 



(1 - vrtd|g„_JJ(gAf,e) 



1 - 



1 - 



De<lM,eDfS\q^jg 



) 



Duqufim J{qM,e) (118) 



It is easy to verify, that if qMfi is a solution to (|^. then AqM,e = 0. 
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Subsequently we calculate the corrections to qm^t using an iterative method for the 
manifold with edge (see Appendix). 

Generally speaking, one could (and should) calculate these corrections for the finite 
models. However, the infinite models are distinguished, because they require such correc- 
tions. 

4.14 The film, and the macroscopic equations 

Let a film be constructed. What next? There are two routes. 

i) Investigation of the conservative dynamics of 'W + 1" variables, where 'W is 
moments, M, and " + 1" is r on the film; 

ii) Derivation of the macroscopic equations for M. 

Actually, the second route is more desirable, it leads to the usual classes of equa- 
tions. The first one, however, is always available, because the film exists always (at least 
formally) but the existence of equations for M is not guaranteed. 

The route of obtaining equations for M is is the same one, suggested by us P5].p^-[T^ 
following Ehrenfest and Zubarev 

i) One chooses a time T. 

ii) For arbitrary Mq one solves the problem of the motion on the film ()65|). ()67|) under 
initial conditions M(0) = Mq, r(0) = tq on the segment t G [0;T]. The solution is 
M(t,Mo). 

iii) For the mapping Mq — > M{T) the system dM/dt = F{M) is constructed. It has 
the property that for its phase flow, 9t{M), the identity 

OriMo) = M{T,Mo) (119) 

is satisfied, for defined T and all Mq. This is a natural projector again (see (jH^ . and 
whole section ). 

In this sequence of actions there are two nontrivial problems: solution to the equations 
on the film, and reconstruction of the vector field by transformation of the phase flow, 
9t, under fixed T. 

The natural method for solving the first problem is the averaging method. The equa- 
tions of motion on the film read 

M = eP{M,T)- f = Q{M,T) (120) 

where e is (formally) small parameter. 

Assuming that the motion of M is slow, one can write down the series of the Bogoliubov- 
Krylov averaging method |SZ] • The first term of this series is a simple averaging over the 
period T : ri(T, M) is solution to the equation f = Q{M, r) under fixed M, 

Mi(t, Mo) =Mo + et (^i £ P(Mo, n{9, Mo))d9^ (121) 
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for t e [0;T], and 

M,{T,Mo)=Mo + e I P{Mo,T{e, Mo))de, (122) 



corresp ondingly. 

The first correction to reconstruction of the vector field, F{M), by the transformation 
of the phase flow, 6x{M), is very simple too: 

F,{M) = ^{eT{M)-M). (123) 
Hence, we obtain the first correction to the macroscopic equations: 

M = F,{M) = ^j^ m{J{qM,r(t,M)))dt, (124) 

where r(t, M ) is a solution to the equation (fHTjl under fixed M (actually, mJ{qM,T) should 
be substituted into (fHTjl instead of M). 

The second and higher approximations are much more cumbersome, but their con- 
struction is not a significant problem. 

Let us demonstrate an explicit expression for M ()124|) for the modeling motion on the 
entropic circle ()51|) . ()112|) for the linear in layer system ()88p . 

The original system is 

/ = JUm{f)) + Lm(f){f - /m(/))' 

where Lm = (/^/J(/))|/;,. 

The macroscopic equations are (see also equations (j53|) ): 

2u} r^"^"^ 

M = m{J{fl,)) + {l/uj)m{LM{q))— sm{ujt)dt 

Jo 

+ (l/cj2)m(LM(gx))— / {I - cos{ujt))dt 

Jo 

= m(J(/;,)) + ^m(LM(g)) + (1/^')(1 - 2/7r)m(LM(g±)), (125) 
where q = Af*^ = J{fli) — 7!'f*^J{fl.j), and Hf*^ is the quasiequilibrium projector (fT^ . 
g = (l-7r^.jLM((l-27r;^)J(/;,)) + DMvr/^m(J(/X,)), 
q± = q- jTTTT — q. 

{\) f*j is the entropic scalar product related with quadratic approximation to the entropy 
fti '■ 

s{f) = s{n,)+Ds\rjf-fij)-ii/2){f-fiAf-fM)n,+o{\\f-n,r), 

{<ii-\qA-)fh 
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Note. In fjl25|) . in accordance with ()67|). the model of motion on the circle ()112|) is taken 
without recalculating the time. Such a recalculation changes the values of the coefficients 
without a change in the structure of the equation: instead of 2/7r and (1 — 2/7r), other 
numbers appear. 

The entropy production for equations (fT^ has the form const '^^^^^c^yature ^^"^^ 

(El. 

In general, equations such as ()125p are determined accurately to the values of the 
coefficients simply by the sequence of the horizon points of the second order finite Kepler 
models, and corresponding gj, g'j 

M = miJiflj)) + Y,{o^MLM{qi)) + A^I^m^))), (126) 

i 

with Oj, (3i > 0. 

The last comment on the positivity of the "kinetic constants" and is important, 
and cannot be easily verified every time. However, in the case under consideration it 
follows from the next theorem. 

The theorem about the positivity of kinetic constants. The motion on the 
Kepler ellipse from start to the horizon point always satisfies the property 

q - qo = aq + /3q; a,/5>0, (127) 

where go is a starting point, q, and q are velocity, and acceleration correspondingly. 

This theorem follows from elementary theorems about analytical geometry of second- 
order curves. 

For the modeling motion on the circle, strictly speaking, this is not so every time. 
Positivity of the coefficients is guaranteed only for m{L{q)), and m{L{q±)). 

Two phenomena can be related to the increase of the number of terms in (jl26p : i) 
alteration of the kinetic constants (terms are not orthogonal to each other, therefore, new 
terms contribute to the previous processes), n) appearance of new processes. 

Motion on an infinite film can lead to the stabilization of kinetic constants as functions 
of M, but it can also lead to their permanent transformation. In the second case one has 
to introduce into macroscopic equations an additional variable, the coordinate, r, on the 
film. 

From the applications point of view, another form of equations of motion on the 
film could be more natural. In these equations kinetic coefficients are used as dynamic 
variables. Essentially, this is just another representation of equations (jU7|) . For every 
kinetic coefficient, k, expression dk/dt = ipki^, M) = ipk{k, M) is calculated in accordance 
with (jnSI), (inZj)- Substitution of variables (r, M) — ^ (fc, M) in this equation is possible (at 
least locally) if value k does not stabilize during the motion on the film. Finally, we have 
the system in the form: 

M = m( J(/;,)) + J2 kjF.iM); k, = ^^{k,, M). (128) 
j 

For the motion starting from the quasiequilibrium state the initial conditions are kj = 0. 
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4.15 New in the separation of the relaxation times 

The classical Bogoliubov's concept about separation of the relaxation times does not agree 
well with the thesis of the quasiequilibrium initial conditions. 

Originally, there are no dissipative possesses in the quasiequilibrium state (the theorem 
of preservation of the type of dynamics for the quasiequilibrium approximation). 

The first thing that occurs during the motion out of the quasiequilibrium initial con- 
ditions is appearance of the dissipation. It can be described (in the first non-vanishing 
approximation) by equation (j^Hj) . It is of special importance here that there is still no 
separation into processes with various kinetic coefficients. This occurs at further relax- 
ation stages: Various processes appear, their kinetic coefficients are determined (see, for 
example, ()126|) ) (or, in some cases, the dynamics of the kinetic coefficients is determined). 
And just after this the "hydrodynamic" relaxation occurs, which is the motion of the 
macroscopic variables to their equilibrium values. 

Generalizing, we can distinguish three stages: 

i) appearance of dissipation; 

ii) branching of dissipation: appearance of various processes; 

iii) macroscopic relaxation. 

It is important to notice in this schema that the determination of the kinetic coeffi- 
cients can occur at both stages: at the second stage when macroscopic (hydrodynamic) 
relaxation can be described in the usual form with kinetic coefficient as functions of the 
macroscopic parameters, as well as in the third phase (motion on the film), when the 
hydrodynamic description includes dynamics of the kinetic coefficients also. 

5 Conclusion 

To solve the problem of irreversibility we have introduced the notion of the macroscopically 
definable ensembles. They are result of evolution of ensembles from the quasiequilibrium 
initial conditions under macroscopic control. The quasiequilibria (ensembles of conditional 
maximum of the entropy under fixed macro-variables) are intensively used in statistical 
mechanics after Jaynes Papers of Rosonoer and Kogan 0S]-|17] significantly affected 
our initial investigation. The primitive macroscopically definable ensembles appear as 
results (for t > 0) of motions which start from the quasiequihbrium state (at t = 0). 
The hypothesis of the primitive macroscopically definable ensembles is very important 
from constructive point of view: Any macroscopically definable ensemble can be approx- 
imated by primitive macroscopically definable ensembles with appropriate accuracy. In 
accordance to this hypothesis it is possible to study the one curve for every value of 
macroscopic variables. These curves form the film of nonequilibrium states. 

The hypothesis about the primitive macroscopically definable ensembles is real hypoth- 
esis, it can be true or false. There is the significant difference between this hypothesis and 
the thesis about macroscopically definable ensembles. The thesis can be accepted, or not. 
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but nobody can prove the definition, even the definition of the macroscopically definable 
ensembles. 

Technically, the solution to the problem of irreversibility looks as follows: we can 
operate only with the macroscopically definable ensembles; the class of these ensembles 
is not invariant with respect to the time inversion. The notion of the macroscopically 
definable ensembles moves the problem of irreversibility into a new setting. It could be 
called a control theory point of view. The key question is: Which parameters can we 
control? It is those parameters that are fixed until "all the rest" come to an equilibrium. 
The quasiequilibrium states are obtained in such a way. 

The further development of this direction must lead to an investigation of the macro- 
dynamics under controlled macro-parameters. This will be a supplement of the postulated 
quasiequilibrium initial conditions with an investigation of a general case of an evolution 
of the controlled ensembles: The initial condition is quasiequilibrium, after which one 
carries out on the system by available control infiuences. 

The method of the natural projector allows us to construct an approximate dynamics 
of macro-variables. Under tendency of the time of projection, r, to infinity, these equa- 
tions should tend to the actual equations of macro-dynamics, if the latter exist. This a 
hypothesis about their existence for the thermodynamic limit (first, the number of par- 
ticles — > oo, and then, the time of projection r — oo) is the basis of the Zubarev 
statistical operator Here, however, we need to make a note. Frequently, physicists use 
objects whose existence and unicity are not proven: solution to the hydro- gaso-dynamics, 
kinetic equations etc. Often, the failure to prove the theorems of existence and unicity is 
treated as an absence of an adequate mathematical statement of the problem (definition 
of spaces etc.). For all this, it is assumed, that substantial obstacles either are absent, or 
can be distinguished separately, independently on the theorem proof in physically trivial 
situations. Existence (or non-existence) of the macroscopic dynamics is a problem of an 
absolutely different kind. This question is substantial: the cases of non-existence can be 
found as frequently as the usual existence. 

The notion of the invariant film of non-equilibrium states, and the method of its ap- 
proximate construction allow us to solve the problem of macro-kinetic even in cases when 
there are no autonomous equations of macro-kinetic. The existence of the film appears 
to be one of the physically trivial problems of existence and uniqueness of solutions. 

Using the Taylor's expansions of the natural projector, the first applications already 
have been constructed [211 113 US]- Particularly, the post-Navier-Stokes hydrodynamic, 
replacing the Burnett equations, have been found. It is free from unphysical singularities 

[2niiii|. 

The formula for entropy production 

defect of invariance 

a = const 

curvature 

makes the geometrical sense of dissipation clear. 

Nevertheless, at least one important problem remains unsolved. This is a problem of 
undivisible events: For macroscopically small time small microscopic subsystems can go 
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through "the whole hfe", from the beginning to the hmit state (or, more accurate, to 
the hmit behaviour which may be not only a state, but a type of motion, etc.). The 
microscopic evolution of the system in a small interval of macroscopic time can not be 
written in the form 

A/ = fAt, 

if it is really the system with microscopic structure, and consists of a large number of 
microscopic subsystems. The evolution of microscopic subsystems in a macroscopically 
small time At should be described as a "ensemble of undivisible events" . An excellent hint 
gives us the Boltzmann equation with undivisible collisions, another good hint gives the 
chemical kinetics with undivisible events of elementary reactions. The useful formalism for 
description such ensembles of undivisible events is developed. It is the "quasi-chemical" 
representation (see elsewhere, for example, the review |221)- But the way from a general 
systems to such ensembles remains unclear. It is a challenge to the following works. 

6 Appendix 

6.1 The method of invariant manifolds 

The aim of this appendix is to give a short presentation of the method of invariant 
manifold, including positive-invariant manifolds with the fixed edge. 

6.2 Construction of the invariant sections 

Let E he a vector space, in the patch U G E the vector field (microscopic system) 

/ = J(/), (/ G U). (129) 

is defined. 

J is assumed to be smoothly continued to the closure of U positively invariantly with 
respect to ()129|) . It means that every solution to ()129|) . f(t), starting under t = in f/, 
is to be found within U for every t G [0; +oo). 

Let S be a vector field (of macroscopic variables), and a surjective mapping m : E ^ B 
is defined. 

It is required to construct such a mapping 

M ^ f* (M G m(f/), f* G f/), (130) 

that Tn{ffj) = M, and is a positive invariant manifold of the system (jl29|) (since U is 
positive- invariant, it is sufficient to verify a local condition: the field J{ffj) is tangent to 
the manifold for each M G m{U)). 

Actually, we continue to keep such a level of strictness (unstrictness) of reasoning when 
such details as topology in E and B etc. are not discussed. If necessary, it could be made 
for particular realizations. 
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We call the mapping ()130|) a section, and the problem of construction of positive 
invariant manifold, an invariant section problem. 

It could be solved with many methods. Here are some of them: 

i) the Taylor expansion over the degrees of an appropriate parameter in a neighborhood 
of the initial approximation (for example the Chapman- Enskog method 

ii) the Newton method (as in the KAM-theory jSU], jlHI but with a incomplete lin- 
earization, as in the original formulation of the method of invariant manifolds for dissi- 
pative systems 

iii) the implementation of Galerkin approximations for each iteration. 

Omitting the well known expansions of the perturbation theory, we consider the direct 
methods. 

In order to make one step of the Newton method with incomplete linearization, we 
need: 

i) an approximate manifold, which we call f2; 

ii) projector P mapping a neighborhood of on Vt. 

For each fT, the projector, 7r,#, mapping E on the tangent space, T.#f2, is 7r.# = 

JM 1 M 1 M 

DfPn\r# (tt is differential of Pq). 

■' Hi 

Usually, the projector is defined in such a way that the layers (prototypes of points 
under projection) could be patches on affine subspaces of E. 
For each we define the invariance defect 

^f*=JUti)-^f*JUti)- (131) 

J M J M 

The invariance equation 

A,#=0 (132) 



is solved with the Newton method with incomplete linearization: for every M we search 
for 5 r# , such as: 

JM 



Pniffi + Sfu) - fti^ 
(l-TT #)DJ(/)| = -A # 

J M J M J M 



(133) 



M- 



If the layers of Pq are patches on the affine manifolds, then ()133|) is a system of linear 
equations. Another form of this system is 



(l-^,*)/^J(/)|,*//* = -A,*. 



(134) 



It should be stressed that in equation ()133|) incomplete (in contrast to the Newton 
method) linearization has been used. We did not differentiate 7r.# in A,# in (I131D . 

Jm Jm 

The discussion is given in (jHJ. We note only that for the simplest stable self-adjoin lin- 
ear systems with incomplete linearization equations (|133|) lead generically to an invariant 
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subspace with the largest (i.e. closest to zero) eigenvalues. In contrast, procedures with 
whole lonearization lead in this case to the subspace closest to the initial approximation. 

As soon as 6f^ is found from equations ()133|) . we substitute for /m+^/m' construct 
new projectors and repeat the procedure. 

Solution to the invariance equation ()132j) by the Newton method with incomplete 
linearisation can turn into a difficult problem. In spite of their linearity, equations ()133|) 
cannot be solved easily every time. One can try to simplify the problem passing from 
the invariance equations to the Galerkin approximations. The simplest example is one- 
dimensional approximations when dfZr = ^M^t* , and equation is solved in the projection 

JM 

on A : 

{A.#\{1 - 71 .#)J{f* + Sm A f*)) f# =0. (135) 

JM JM JM JM 

The entropic scalar product (1) .# = —D^S(f )\r# is used. 

Solving (jl35p with the Newton method, we obtain at the first iteration: 

(A ,# I A ,# ) r# 

J)M = /a/ Jm Jm /ioc\ 

{A.#\DJif).#A.#).#- ^ ^ 

JM JM JM JM 

More often for the dissipative systems the denominator is negative, and this allows us 
to move on. For the conservative systems the one-dimensional Galerkin approximations 
lead to an unsatisfactory result, at least, in the combination with the Newton method 
(with the incomplete linearization). 



6.3 The entropic thermodynamic projectors 

The simplest choice of Pn is obvious: 

Pnif) = ft: (137) 

for each value of / the values of macroscopic variables, m{f), can be calculated. Based 
on these values, the corresponding point with the same value of m(/) on the manifold 
can be obtained. 

However, projector (jl37|) does not satisfy the physical constrains every time. 

On the set U a concave function, the entropy, S, is defined. Two kinds of systems 
are under consideration: i) (jl29j) dissipative, for which dS/dt < 0, in accordance with the 
system, ii) conservative, for which dS/dt = 0. The quasiequilibrium manifolds, are 
physically distinguished. They are the solution to the variational problem 

S{f) max, 

m(/) = M. (138) 

Application of the simplest projector ()137p leads to the fact that the vector field, tt/.^ J(/^), 
preserves the type of dynamics of the system on the quasiequilibrium manifold. For the 
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conservative J it is conservative, and for the dissipative systems it is dissipative too (with 
the same entropy). Such a preservation of the type of dynamics by the projector is 
guaranteed only for the quasiequihbrium manifolds. 

However, practically for every manifold Vt = {/^} it is possible to construct such a 
projector, Pq, that every is a solution to the problem 



in a neighborhood of f2. For this it would be sufficient that for every M the functional 



Thus, in order to control the physical sense of obtained approximations, one needs to 
consider the projector depending on the manifold PT] . 

6.4 Method of invariant manifold for the positively invariant 
manifolds with fixed edge 

In the problem of construction of invariant sections the position of points was not 
fixed. Only fulfillment of condition m{ffj) = M was important. There is another kind 
of problems where one needs to find a positive invariant manifold with a fixed edge. 
Practically, this is a problem of construction of a trajectory of the edge in accordance 
with (fT^ for t e [0;+oo). 

These problems include the problem of initial layer ^H]; as well as the problem 
of construction of the film of non-equilibrium states discussed in this paper. 

The iterative methods described above cannot be implemented here, because they 
destroy the boundary conditions on the edge of the manifold. If the invariance conditions 
are fulfilled on the edge of the initial approximation, f^, accurate to the A;— th derivative 
in time, then the Newton method leads to the fact that after the {k + 1)— th iteration the 
edges of the manifold will be changed (the film tears off the edge) . 

In the previous paper devoted to the problem of initial layer [15 we have technically 
overcame this difficulty. To do this, we simulated the trajectory as an elastic beam with a 
rigidly fixed end. In the mechanical equilibrium this beam had the form of an approximate 
trajectory. Later, it was elastically attracted to the result of the Newton iterations. Even 
though that this technique allows to avoid the separation of contact between the edge and 
the film, its technical artificiality forces us to continue to search for new methods. 

The application of the Picard iterative procedure allows to conserve initial conditions. 
For the film equation we write: let q^^ ^ be an approximation for the film, then the Picard 
iteration gives: 




(139) 



DS(f)\f# eliminate ker 7r.# : 

Jm 1m 




(140) 




(141) 
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Here dqM,@/d6gO^ ^ is the right hand side of (j3Hj) taken at the point q%ifi. 

Let us define, as usual, /S.m,b = dqM,e/d9\qO^^ —dq% .^/d9 as the difference of the vector 
field and its projection on the approximate manifold, then the Picard iteration obtains 
the form: 



gl,,, = qlj^, + / AM,ede. (142) 
Jo 

The Picard iteration gives a good result for small r but can be too radical for large 
one. It is possible to use the Picard iterations together with the weight which ensures 
essential dependence of the correction, g^/^, not from all AM,e, < 9 < t, but from those 
within a segment 6 E [r — h^r]. For example, 

QM,r = <lM,r + f C^' ' ^ A m fidO . (143) 

Jo 

Another choice of the weight function is possible. For large r and sufficiently small h 
and A 7^ 0, the formula ()143|1 gives: 



d/,r = &r + ^AAf,r + 0(/l). (144) 

For small r we obtain: 

= qll,r + -^^M^r + o{t A) . (145) 

where k is order of zero of A^^^ at the point r = 0. 
Joining ()144|1 and ()145j) . we obtain: 

In all cases the question of how to choose the step arises. The simplest solution exists 
for ()146|1 : it is possible to take the step, h, depending on the point: 

Qll,r ~ QM,r + ^M,rAM,r. (147) 

where Am,t = min{r/(A; + l), 5m,t}^ and 5m,t is to be found from the condition of stopping 
in the direction of A ()135|) . 

Various combinations of the Picard and Newton iterations can generate a separate 
subject for investigation. Their simplest hybridisation consists of the following. Let for 
each (M, r) the step, 5m,t, is found according to the Newton method ()133p . 

Define 

Qkr = Qli,r + ^/h f e^'-^y'^6qM,ed9. (148) 

In order to have the step value on the direction 6qM,e close to 1 for large r, the 
multiplier 1/h has been used. An analog to (jl46p is 
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The typical time scale, which separates in ()137p and ()138p the Picard (small times, 
T <^ h), and Newton (large time, t ^ h) iterations, can be estimated by the curvature 
radius (in general, from the relations ||g||/||g|| or llfjU/UgxH ). 

Further development of the methods will be determined by particular applications. 
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